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ABSTRACT 


We present results on the dust attenuation curve of z ~ 2 galaxies using early observations from the 
MOSFIRE Deep Evolution Field (MOSDEF) survey. Our sample consists of 224 star-forming galaxies 
with nebular spectroscopic redshifts in the range Zgpec = 1-36 — 2.59 and high S/N measurements of, 
or upper limits on, the Ha and H/3 emission lines obtained with the MOSFIRE spectrograph on 
the Keck I telescope. Using deep multi-wavelength photometry, we construct composite spectral 
energy distributions (SEDs) of galaxies in bins of specific star-formation rate (SFR/M*) and Balmer 
optical depth. These composites are used to directly constrain the shape and normalization of the dust 
attenuation curve over the full wavelength range from the UV through near-IR for typical star-forming 
galaxies at high redshift (z > 1.4). Our results imply an attenuation curve that is very similar in shape 


and normalization to the SMC extinction curve at wavelengths A > 2 500 A. At shorter wavelengths, 
the shape of the curve is identical to that of the Calzetti et al.j (20001 starburst attenuation relation, 
but with a lower normalization (Rv), implying less attenuation at a fixed wavelength for a given SED 
shape. Hence, the new attenuation curve results in SFRs that are « 20% lower, and stellar masses 
that are Alog{M*/M q) ~ O.lGdex lower, than those obtained with the starburst attenuation curve. 


We find marginal evidence for excess absorption at 2175 A. Moreover, we find that the difference in the 
reddening—and the total attenuation—of the ionized gas and stellar con tinuum correlate s strongly 


with SFR, such that for dust-corrected SFRs > 2OM0yr assuming a Chabrier (2003) IMF, the 


nebular emission lines suffer an increasing degree of obscuration relative to the continuum. A simple 
model that can account for these trends is one in which the UV through optical stellar continuum is 
dominated by a population of less reddened stars, while the nebular line and bolometric luminosities 
become increasingly dominated by dustier stellar populations for galaxies with large SFRs, as a result 
of the increased dust enrichment that accompanies such galaxies. Consequently, UV- and SED-based 
SFRs may underestimate the total SFR at even modest levels of « 20 Mq yr~^. 

Subject headings: dust, extinction — galaxies: evolution — galaxies: formation — galaxies: high- 
redshift — galaxies: star formation 


1. INTRODUCTION 

As a fundamental byproduct of star formation, dust 
plays an important role in our understanding of heavy 
element synthesis in stars, the dispersal of such elements 
into the ISM, and our interpretation of the stellar popula¬ 
tions of galaxies given the wavelength-dependent effects 
of dust absorption. The latter aspect is encoded in the 
shape of the attenuation curve (jDraine & Li|2007), which 
conveniently parameterizes the complexities of the dust 
properties and spatial distribution of that dust with re- 
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spect to the stars in galaxies. Observations of the Milky 
Way and nearby galaxies have yielded detailed informa¬ 
tion on the shape of the dust extinction and atte nuation 
curves, respectively, in the local Universe (e.g., Prevot 
et al.||19^ Cardelli et al. 1989 Pei||1992 (Jalzetti et al. 


19941 120001 lOordon 

WilJ et al.lp!^ 


et al.l 2003 Johnson et al. 2007a 

Armed with such information, one 
can then deduce the amount of attenuation at any given 
wavelength (i.e., the dust obscuration) and recover the 
bolometric star-format ion rate (SFR), accounting for the 
starlight absorbed and re-radiated by dust. 

Understanding the shape of this curve takes on added 
importance for evaluating the global energetics and stel¬ 
lar populations of galaxies at high redshift, where di¬ 
rect measures of dust emission from far-infrared (far-IR) 
data become observati onally inaccessible f or typical star¬ 
forming galaxies (e.g., Reddy et al.||2006 ). The stacking 
of mid- and far-IR data from the Spitzer Space Telescope 
and the Herschel Space Observatory has proved benefi¬ 
cial for our understanding of the average du st properties 
of typical L* star-forming galaxies at z ~ 2 ( Reddy et al. 
2006| |Daddi et al.||2007t |Reddy et al.||2012a!|). IR obser¬ 
vations of lensed star-forming galaxies can lend insight 
into the dust emission from individual L* (and fainter ) 
galaxies (e.g., Siana et al.|2008 2009 Sklias et al.|2014 ). 
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but these samples remain small. These stacking and 
lensed galaxy studies lack the precision required to un¬ 
cover changes in dust properties and dust/stars geometry 
from galaxy to galaxy. Further, while the UV contin- 
uum slope (e.g., Calzetti et al.||1994 Meurer et al.||1999 l 


has been shown by numerous studies to be a reasonable 
proxy for dust attenuation in star-for ming galaxies in the 
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tor individual galaxies to very faint magnitudes at 


high redshift ( 
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2012 Oesch et al. 
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IReddy et al.||l: 

iolOl lUo- 

quien et ai.|2012| Wilkins et al. 

|2013 see also o 

elow). In 
be com- 

addition, the measurement of t 

he U V slope car 


|Noll et al. 

2009 

Buat et al. 12011 

Kriek & Gonroy|2013 

Buat et al. 

2012 

. Despite these issues, the U V slope is 


shift given that it can be measured easily for individual 
star-forming galaxies, particularly those selected by their 
rest-frame UV colors. 

Separately, existing studies of the shape of the at¬ 
tenuation curve at z > 1.5 have relied either on spec¬ 
troscopy over limited wavelength ranges, indirect infer¬ 
ences based on SED modeling of br oad- and medium- 


band data and/or SFR comparisons i 

Reddy et al.| 
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Erb et ai.|2006al|Noll et al. 120091 IRedc 

ly et al.|2010l|Buat| 

et al. 201H |2012 Kriek ik, Conroy 2013P, or have been 


tral shape (e.g., galaxies with weak Balmer and 4000 A 
breaks) and UV luminosities with an explicit assump¬ 


tion for th e shape of the intrinsic spectrum (Scoville 
et al. 20151. What have been needed for a more robust 


study of the attenuation curve are indicators of dust and 
SFR at high redshift that: (a) can be measured for in¬ 
dividual L* galaxies at high redshift; (b) are less atten¬ 
uated than the UV luminosity; (c) are less sensitive to 
age and star-formation history; and (c) probe star for¬ 
mation on the shortest timescales. When coupled with 
samples of galaxies with near-IR spectra and broad- and 
intermediate-band photometry spanning the UV through 
near-IR wavelengths, these SFR/dust indicators can be 
used to address the shape of the attenuation curve of 
high-redshift galaxies. 

The HI recombination line luminosities and their ratios 
have been considered the “gold standard” for m easur¬ 


ing SFRs and dust a ttenuation in local galaxies (Ken- 
nicutt & Evans 2012). In particular, the Balmer (and 


higher series) transitions in star-forming galaxies arise 
from the HII regions generated by the most massive O 
stars (> 8Mq) and therefore trace the SFR averaged 
over timescales of just a few Myr, making these trac¬ 
ers insensitive to the star-formation history on longer 
timescales. In nearby galaxies, these lines are easily ob¬ 
servable from the ground to faint limits, and they are 
less affected by attenuation than UV emission. The 
higher series lines, including Paa (A = 1.875 /im) and 


Bry (A = 2.166 ^m), are particularly useful indicators 
as their relatively long wavelengths imply negligible dust 
absorption. However, these features are weak relative 
to the Balmer lines and, at z > 1.5, shift into the mid- 
IR where current instrumentation has been insufficient 
to detect the lines for all but the most luminous and 


strongly gravitationally lensed galaxies (jPapovich et al. 

200 ^. 


Until recently, obtaining measurements of nebular re¬ 
combination (other than Lya) has been prohibitive for 
faint galaxies at redshifts z > 0.5 as these lines are 
shifted into the near-infrared (near-IR) where the terres¬ 
trial background is substantially higher than in the opti- 


narrowband imaging (IGeach et al. 

20081 

Lee et al.||2012 

Ly et al.|2012 

Sobral et al.|2012 i' 

domch 

eva et al.|2013). 


hie Space Telescope ( 

HST) 

WFG3 camera (lAtek et al. 

2010| 

Brammer et al.|2012| 

IDommguez et al.|20131|Price 

et al. 

2014p have begun to 

advance the study of nebular 


redshift galaxies (0.5 ^ z < 2.5). Nonetheless, such 
studies are generally limited to detecting high equivalent 
width lines, do not allow for the separation of lines close 
in wavelength (e.g., Ha and [Nil]), or have a low resolu¬ 
tion (i? < 130) with a limited wavelength range that hin¬ 
ders the determination of physically important line ratios 
and widths. Moderate-resolution near-IR spectros copy 


exists for both UV-selected star-f orming galaxies (Erb 
et al.| 2006b Mannucci et al.|[2009| ) and near -IR-selected 


samp les of massive (M > 10^^ Mq) galaxies (Kriek et al. 


2008), yet these samples are quite small ahd typically 


do not cover both the Ha an d H^ lines. In cas es where 
both lines were targeted (e.g., Kriek et al.|2008), the Mj3 


line is commonly undetected given the depths of previous 
spectroscopic campaigns. The advances enabled with the 
new generation of multi-object near-IR spectrographs on 
8-lOm-class telescopes have begun to remedy these past 
deficiencies in high-redshift near-IR sp ectroscopy (e.g., 
Kashino et al.|j2013 Steidel et al.||2014). 

'id this end, we have begun th e ~MOSFIRE Deep 


Evolution Field (MOSDEF) survey (Kriek et al. 2014) 
conducted with the recently-com missioned multi-obj ect 
near-IR spectrograph MOSFIRE (McLean et al. 2012) on 
the Keck I telescope. MOSDEF aims to characterize the 
SFRs, dust attenuation, chemical content, ISM physical 
conditions, dynamics, stellar populations, and black hole 
accretion activity of a large sample of ^ 1500 H-selected 
galaxies and AGN in the redshift range z = 1.4 —3.8, cor- 
responding to the peak of star formation (see review of 
Madau & Dickinson 201 ^ and active galact ic nucleus 


Silverman et al. 2008). In this 
e of 224 galaxies with measure- 


(AGN) activity (e.g. 
study, we use a samp, 
ments of the Ha and H/3 emission lines from the Hrst 
two years of the MOSDEF survey to provide a direct 
measurement of the shape of the attenuation curve over 
the full wavelength range from the UV through near-IR 
for typical star-forming galaxies at z > 1.4. Our analy¬ 
sis marks a signihcant improvement over existing stud¬ 
ies because: (a) we consider galaxies with a large range 
in spectral shapes (from very blue and dust-free star¬ 
forming galaxies to very dusty and highly star-forming 
ones); (b) we use an indicator of dust, namely the Balmer 
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decrement, which is sensitive to the reddening of the ion¬ 
ized gas and thus is physically distinct from the shape of 
the SED (and thus less sensitive to the details of the star- 
formation history and age of the stellar population), and 
which can be measured for individual galaxies in our sam¬ 
ple; (c) our method of computing the attenuation curve 
is model-independent in that we make no assumption of 
the shape of the instrinic SED—rather, we use the data 
themselves to establish the how galaxy spectral shapes 
vary with increased dust attenuation; and (d) the exten¬ 
sive multi-wavelength photometry affords us the ability 
to measure the shape of the attenuation curve at long 
wavelengths and, thus, the normalization of the curve. 

In Section we introduce the MOSDEF survey and 
the basic measurements including line fluxes, stellar pop¬ 
ulation parameters, and SFRs. Section [^presents a dis¬ 
cussion of the effects of dust on the stellar continuum 
and the derivation of the attenuation curve. The impli¬ 
cations for the attenuation curve on our inferences of the 
stellar populations and SFRs, as well as the total attenu¬ 
ation of the stellar continuum and ionized gas, in distant 
galaxie s are discussed i n Section]^ Throughout we as¬ 


sume a 


Chabrier (2003) initial mass function (IMF) and 
a cosmology with Hg = 70kms“^ Mpc“^, Da = 0.7, and 
Dm = 0.3. Line wavelengths are in vacu um. All mag¬ 


nitud es are expressed in the AB system (Oke & Gunn 


2. DATA AND BASIC MEASUREMENTS 

2.1. MOSDEF Survey 
The MOSDEF survey i s being conducted u s ing the 


MOSFIRE spectrograph (McLean et al. 2012 Steidel 
et al.|[2M4| on the 10 m Keck I telescope. This recently 


commissioned instrument enables the simultaneous near- 
IR spectroscopy of up to « 46 (with typically 30) galax¬ 
ies across a 6' x 3' field of view. MOSFIRE provides a 
wavelength coverage from 0.97 to 2.40 /xm with a spectral 
resolution of i? « 3300, 3650, and 3600 in the J, H, and K 
bands, respectively. Details of the s urvey, observations , 
and data reduction are provided in Kriek et al. (2014), 
and below we summarize the salient aspects relevant tor 
this analysis, which uses the current data from the first 
two observing seasons of the survey. 


2.2. Galaxy Selection 

Catalogs including ground- and space-based photom¬ 
etry from 0.3 to 8.0 /xm, provided to us by the 3D-HST 
grism survey team (PI: van Dokkum; Brammer et al. 
20T^ iSkelton et al.|[20l4| ), were used for target selection 
( [Kriek et al.||2U14 ). 'These catalogs are based on sources 
detected in a weighted combination of the HST F125W, 
F140W, and F160W images in the five fields of the Cos¬ 
mic Assembly Near- Infrared Deep Extragalactic Legacy 
Survey (CANDELS;|Koekemoer et al.|20il||Grogin et al. 
Mlj): COSMOS, ECS, COODS-lNorth, C(50DS-South, 
and UDS. The targeting catalogs include photomet- 
ric redshifts and, where available, HST gris m redshifts 
(Brammer et al. 2012 Skelton et al. 2014) and exter- 
nal spectroscopic redshifts. We restricted the targeting 
catalogs to include only those galaxies with photometric 
and/or grism or external spectroscopic redshifts in the 
three ranges 2 = 1.37— 1.70, 2.09 — 2.61, and 2.95 — 3.80, 
where the various strong nebular emission lines of inter¬ 


est, including [OII]AA3727,3730 A, [OIII]AA4364,5008 A, 
H/3, Ha, and [Nil] AA6550,6585 A, fall in the Y, J, H, and 
K atmospheric windows. Targets are selected to fixed H- 
band magnitude limits of mpieow = 24.0, 24.5, and 25.0, 
for the z = 1.37 - 1.70, 2.09 - 2.61, and 2.95 - 3.80 red- 
shift ranges, respectively, to ensure a roughly consistent 
stellar mass limit of « 10® Mq across the three redshift 
ranges. For brevity, we refer to these as the z ~ 1.5, 

2.3, and 3.3 samples. In this analysis, we consider only 
galaxies from the two lower redshift ranges, where both 
the Ha and H/3 lines are covered. 


2.3. Observations and Efficiency 

Targets were prioritized for spectroscopic observations 
based on two primary criteria. First, bright galaxies were 
up-weighted relative to fainter galaxies to ensure a sta¬ 
tistical representation of rarer, massive galaxies. Sec¬ 
ond, within a given magnitude range, objects were pri¬ 
oritized according to whether they had external spectro¬ 
scopic redshifts, grism redshifts, photometric redshifts 
toward the center of the redshift ranges specified above, 
and photometric redshifts at the edges of the redshift 
ranges, in that order. The sample characteristics, suc¬ 
cess rate, and comparison of the spe ctroscopic sample t o 
the parent sample are presented in Kriek et al. (2014). 
Masks were observed using an ABA'B' dither pattern 
with 0'.'7 width slits for 2hrs each in the J, H, and K 
filters for the z ^ 2.3 sample; and Ihr each in the T, J, 
and H filters for the z ^ 1.5 sample. Our observations 
were taken with a typical seeing of O'.'d — I'.'O. Data re¬ 
duction and spectral extraction were accomplished using 


a custom I DL pipeline and interactive software (Kriek 
et al.||2014). The prioritization and exposure time strat- 
egy have proved highly successful, with a ^ 80% spec¬ 
troscopic confirmation rate on average per mask, and a 
> 70% success rate for objects at the faint limit of our 
z ~ 2.3 sample. Figure shows the spectroscopic red¬ 
shift distribution of 259 galaxies (excluding X-ray and 
IR-selected AGN; Coil et aL||2015 ) in the two lower red¬ 
shift ranges targeted with MUSDFF thus far. 


2.4. Slit Loss Corrections 

Slit loss corrections have a particular bearing on our 
analysis, as they are required to ensure accurate flux ra¬ 
tios between lines taken through different photometric 
filters and observed under different weather conditions. 
Response corrections of the spectra were performed us¬ 
ing telluric standard star observations, and an initial flux 
calibration is implemented by forcing the spectra of a 
bright “slit star” placed on each of the science masks to 
match the star’s broadband photometry. However, the 
typical seeing of our observations is such that most of 
the targeted galaxies are spatially resolved and thus can¬ 
not be treated as point sources when correcting for the 
loss of flux outside the slit apertures. Thus, we modeled 
the F160W light profiles of each of our galaxies using a 
two-dimensional rotated elliptical Gaussian which, when 
convolved with the seeing disk and integrated through 
the slit, yields an estimate of the fraction of flux lost 
outside the slit, relative to the fraction lost for the slit 
star. The efficacy of this procedure is evaluated by com¬ 
paring the spectroscopic magnitudes of galaxies detected 
in the continuum with their broadband magnitudes. Our 
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Figure 1. Left: Spectroscopic redshift distribution of 259 galaxies in the two lower redshift ranges (2 = 1.37 — 1.70 and 2.09 — 2.61) 
targeted with MOSDEF during the first two observing seasons. Right: Ratio of Ha to H/3 emission flux versus observed Ha luminosity, 
for our sample of 259 galaxies with spectroscopic redshifts 1.35 < 2apec < 2.59 and S/Nna > 3. The emission line fluxes are corrected for 
underlying Balmer absorption. The 46 galaxies with 3cr upper limits i n H/3 flux are den oted by the triangles. The horizontal blue line 
indicates the intrinsic ratio (Ha/H/3 = 2.86) for typical ISM conditions (|Osterbrock|1989||. Of the 213 galaxies with significant detections 
of both Ho and H/3, 58 and 85 are consistent with having non-zero nebular attenuation at the 3 cr and 2 cr levels, respectively. 


slit loss correction procedure yields an absolute spectral 
flux calibration that is typically within 18% of that pre¬ 
dicted from the photometry, with a dispersion of roughly 
the same order. 

2.5. Line Flux Measurements 

Line measurements were performed by first estimat¬ 
ing an initial redshift based on the highest signal- 
to-noise (S/N) line for each galaxy, and then fitting 
Gaussian functions at the observed wavelengths of all 
other lines. We allowed for a linear fit to the con¬ 
tinuum underlying each line, and the [OIIJA3727,3730 
and Ha-|-[NII]A6550, 6585 lines were fit using double and 
triple Gaussian functions, respectively. To estimate line 
flux errors, the spectra of each galaxy were perturbed 
1000 times according to their error spectra and all the 
lines were remeasured from these realizations. 

The Ha and H/3 line fluxes were corrected for under¬ 
lying Balmer absorption arising primarily from the at¬ 
mospheres of A stars. As the continuum is generally 
undetected in the MOSDEF spectra, we estimated the 
Balmer absorption from the stellar population model fits 
to the galaxies, as discussed in the next section. For each 
model fit, we measured the absorbed flux relative to the 
continuum. Generally, the corrections reach peak rest- 
frame equivalent widths of < 3 A and < 4 A for Ha and 
H/3, respectively, at an age of 500 Myr. For a constant 
star-formation history, the equivalent widths equilibrate 
to « 2 — 3 A for older ages. The typical corrections to 
the Ha and H/3 emission line fluxes are < 3% and < 10%, 
respectively. The corrected Ha/H/3 ratios for our sample 
are shown in Figure [l] The random uncertainties in the 
Balmer absorption line fluxes were estimated by calculat¬ 
ing the dispersion in the absorbed fluxes for model stellar 
population fits to the photometry perturbed many times 
according to the photometric errors. The uncertainties 
in the absorbed fluxes are ~ 2% and much smaller than 


the typical « 15% uncertainties in the emission line flux 
measurements, and therefore we chose not to include the 
additional random error of the absorbed fluxes in our 
analysis. 

It is useful to assess whether the 46 galaxies with 3 a 
upper limits in H/3 have distinctive properties relative 
to those of the 213 H/?-detected galaxies. The 46 H/3- 
undetected galaxies have lower limits in Ha/H/3 that 
span the full range observed for the H/l-detected galax¬ 
ies, and the broadband SEDs for the former are otherwise 
unremarkable compared to the latter. Moreover, the fre¬ 
quency with which the H,0 line is significantly affected 
by strong skylines is much higher (37 out of 46, or 80%) 
for the H/3-undetected galaxies than for the H,5-detected 
galaxies (60 out of 213, or 28%). Thus, the H,5 non¬ 
detections can at least partly be attributed to the fact 
that the galaxies lie at redshifts where H/3 happened to 
fall on or close to a strong skyline. 

To further investigate any systematic differences be¬ 
tween these 46 galaxies and the rest of the sample, we 
constructed weighted average spectra of the H/3-detected 
and -undetected galaxies and measured the Ha and H/3 
fluxes in these combined spectra. The weighting of in¬ 
dividual spectra when computing the average is simply 
1 /ct(A)^, where ct(A) if the value of the error spectrum at 
wavelength A. The H/3-detected and -undetected galax¬ 
ies have (Ha/H/3) =4.1 and 4.3, respectively, after cor¬ 
recting for the mean Balmer absorption in these galaxies. 
This difference is smaller than the dispersion in the mean 
Balmer decrements obtained when we bootstrap resam¬ 
ple the galaxies contributing to the mean stacks. The 
galaxies with undetected H/3 have an average Balmer 
decrement that is similar to that of the H,5-detected 
galaxies, implying that they suffer the same degree of 
obscuration on average. Only those galaxies with > 3cr 
detection of H/3 (and Ha) are considered in the derivation 
of the attenuation curve, but we specifically include the 
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sample of H/3-undetected galaxies throughout the paper 
as necessary. 


2.6. Stellar Population Modeling 

To aid our analysis, we modeled the stellar populations 
of galaxies in our sample to derive SED-based SFRs, 
color excesses, ages, and stellar masses. The basic proce¬ 
dure was to model the multi -wavelength photometry of 
the galaxies using a suite of Bruzual & Chariot (20031 
stellar population models for ditterent star-tormation 
histories, ages, and r eddening; the latter assumes the 
Calzetti et al. (2000) attenuation curve, but we con- 
sider other attenuation curves below. The photometry 
is corrected, if necessary, for the contribution from the 
strongest emission lines measured in our survey, includ¬ 
ing Ha, H/3, [OIII]AA4960,5008, and [OII]AA3727,3730. 
We modeled the stellar populations using constant, de¬ 
clining, and exponentially rising star-formation histories 
with e-folding times of r = 100 — 5000 Myr, and gen¬ 
erally adopted the exponentially rising histories as the 
default for reporting results from the fitting. In the 
case of 14 galaxies, exponentially declining models were 
adopted as they yielded better fits to the rest-UV pho¬ 
tometry. Incorporating the parameters obtained for a 
constant star-format ion history does not alter our anal¬ 
ysis or conclusions. In the modeling, we allowed the 
age to vary fro m 50 Myr (approxima tely the dynamical 
timescale; e.g., Reddy et al. 2012b I to the age of the 
Universe at the redshift of each galaxy. Reddening in the 
range 0.0 < E{B — V) < 0.6 was considered. We adopted 
the parameters corresponding to the model that gave a 
minimum relative to the photometry. Errors in the 
best-fit parameters were derived by perturbing the pho¬ 
tometry assuming a Gaussian distribution of full-width- 
half-max equal to the photometric uncertainties. Each 
of these realizations of the photometry were then refit to 
find the best-fit model, and the errors in best-fit param¬ 
eters were taken to be the values that encompassed 68% 
of the simulations. 

Of the 259 galaxies in our sample, 35 (31 H/3-detected 
and 4 H/3-undetected) galaxies had noisy photometry 
relative to the best-fit model and were excluded from 
subsquent analysis. Thus, the final sample consists of 
224 galaxies, with 182 H/3-detected galaxies and 42 H/3- 
undetected galaxies. 


2.7. Balmer Optical Depth 


As in Calzetti et al. (1994), we define the Balmer op¬ 
tical depUi 


V 2.86 ) 


( 1 ) 


2.8. SFR Calculations 


The default measurements we use for the SFRs come 
not from the SED fitting, but from the Ha flux mea¬ 
surements. The latter are more immune to the degen¬ 
eracies inherent in SED modeling, less affected by ex¬ 
tinction than the r est-UV, and trace the “instantaneous” 
SFR. We used the Kennicutt (1998) re lation to convert 
Ha luminosities to Sl'Rs, assuming a Chabrier (2003) 
IMF. These observed SFRs were then corrected tor at- 
tenuation by assuming an attenuation curve (fc(A)) and 
nebular color excess {E{B — U))gas)- In most previous 
investigations at high redshift, the attenuation curve as¬ 


sumed for the nebular line regions is the Calzetti et al. 
(2000) curve, with color excesses based o n those found 


tor th e stellar continuum (also assuming a Calzetti et al. 
2000 curve). In our analysi s , we calculated ShK(Ha 


assuming the 
curve, as was 


Cardelli et al. (1989 ) Galactic extinction 
he original intent of 


alzetti et al. (1994), 
but we explore the effect of assuming other attenua- 
tion curves in deriving the color excess, as well as the 
implications for dust corrections and bolometric SFRs, 
in Sections |4.3| and 14.4| below. Throughout the paper, 
and unless noted otherwise, the specific SFRs (sSFRs) 


are computed assuming the dust-corrected SFRs(Ha) 
and the SED-determined stellar masses. This approach 
is more advantageous than combining SED-determined 
SFRs with M* as these quantities are highly correlated 
with each other (i.e., both SFR(SED) and M* are con¬ 
strained by the normalization of the best-fit SED model 
with the photometry). 


3. THE EFFECT OF DUST ON THE STELLAR 
CONTINUUM EMISSION 

As noted above, the distribution of Ha/H/3 ratios im¬ 
plies the presence of dust in the majority of the objects in 
our sample. Before considering the level of obscuration 
required to reproduce the observed Balmer decrements— 
a step that requires the assumption of an attenuation 
curve—it is instructive to examine how the galaxy spec¬ 
tral shapes vary with Ha/H/3 ratio. We focus specifically 
on the sensitivity of the UV spectral slope to dust obscu¬ 
ration, given the large absorption cross-section of dust to 
UV photons, and the prevalent use of the UV slope as 
a dust indicator for typical star-forming galaxies at high 
redshifts. Formally, the UV slope /3 is defined such that 
f\ oc A^. Note that /3 is the observed UV continuum 
slope, irrespective of the underlying stellar population. 
Thus, a galaxy with an older stellar population and a 
significant contribution of UV flux from stars of spectral 
type A and later, can exhibit the same spectral slope as 
a galaxy whose UV continuum is dominated by OB stars 
and contains significant amounts of dust. 


Here, Tb is the difference in optical depths measured at 
the wavelengths of H/3 and Ha, with the usual defini¬ 
tion: T\ = — where I\ and are the observed 

and intrinsic intensities, respectively. Thus, Tb is directly 
related to the color-excess, E{B — V). Objects with mea¬ 
sured Ha/H/3 < 2.86, the theoretical minimum value 
in the absenc e of dust for Ca se B recombination and 
T = 10000K ( |Osterbrock||l989| ), are assigned Tb = 0|^ 


^ The intrinsic ratio Ha/H/3 is insensitive to the electron density, 
varying by less than 2% from rie = 10^ to 10® cm“®. 


3.1. Calculation of the UV Slope 

In practice, /3 can be measured directly from a high 
signal-to-noise spectrum using a number of continuum 
windows spanning the range 1250 — 2600 A, where the 
windows are designed to exclude regions of strong inter¬ 
stellar/stellar absor£tian_Md_t^ 2175 A dust absorp¬ 
tion feature (|Calzetti et al. 1994). For most galaxies 
lacking rest-UV spectroscopic data, as in our case, ei¬ 
ther broad- and intermediate-band photometry, or the 
stellar population model that best fits such photometry. 
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can be used to estimate /3. We calculated a “photomet¬ 
ric” UV slope (/3phot) by fitting directly the photometry 
(taking into account the photometric errors), which in¬ 
cludes for all objects in our sample anywhere from 4 to 
19 photometric points that lie in the (rest-frame) range 
A = 1250 — 2600 A. The measurement uncertainty in the 
rest-UV photometry results in best-fit UV slopes with a 
typical error of A/3phot — 0.2 
For comparison, we also computed an “SED” UV slope 
(/^sed) by fitting the flux of the best-fit stellar population 
model as a function of wavelength, using only those flux 
p oints that lie w i thin t he 10 continuum windows given 


Calzetti et al. (1994). The typical formal uncertainty 


in /tsED, when using the 10 aforementioned windows, is 
A/1sed — 0.1 (not including the systematic error associ¬ 
ated with the assumed stellar population model). The 
two measures of UV slope, /3phot and /3 seDi are highly 
correlated because the stellar population model is fit to 
the same photometry used to calculate /3phot- 

3.2. Relation between UV slope and Balmer Line 
Opacity 

Our expectation is that if the UV slope is sensitive 
to dust obscuration, which a number of studies have 
shown to be the case for different subsets of high-redshift 


Reddv et ai.|2UU6l 

Uaddi et al. 

5I1U7 

Fannella et al. 


Buat et ai.|2009| Magdis et al. 

2010 

Reddy et al.|2012ap 


then we should observe a correlation between /3 and Tb. 
The scatter between these two quantities may hint at 
other factors that influence /3, including star-formation 
history, metallicity of the underlying stellar population, 
and differences in the geometrical distribution of stars 
and dust in galaxies. 

Figure shows the distribution of /3phot with Tb for 
the 182 Ejects in our sample with S/N's_p > 3. The 
average Ha/H/3 ratio for the 42 H/3-undetected galaxies, 
as calculated from the composite spectrum (see above), 
is 4.1 ± 0.1, corresponding to Tb = 0.37 ± 0.02. The 
average value of the individual /3phot measurements for 
these 42 galaxies is /3phot = —1.00. These average val¬ 
ues for the H/?-undetected galaxies are indicated by the 
star in Figure [^. Generally, there is a large dispersion 
in /3phot at a given Balmer optical depth. In the limit 
where Tb —>■ 0, /3 varies from the bluest (/3phot — —2.3) to 
reddest (^Sphot — —0.5) values observed. As noted above, 
P can be sensitive to factors other than the dust con¬ 
tent, possibly accounting for some of the scatter in UV 
slope at a given Balmer optical depth. These factors are 
discussed in detail below. 


3.3. Effect of Star-formation History on the UV Slope 
3.3.1. Variations with Specific Star-Formation Rate 

Galaxies with significant older stellar populations can 
exhibit UV slopes that are redder than those predicted 
from dust attenuation alone. Specifically, A stars from 
previous generations of star formation can contribute sig¬ 
nificantly to the near-UV flux if the current SFR—and 
hence the contribution from O and B stars—is low, as 
might be the case for galaxies with low specific SFRs 
(sSFRs). The two consequences of this effect are to in¬ 
crease the scatter in /3 at a fixed Tb, and to accentuate 
the trend in /? vs. Tb if dustier galaxies are on average 


Table 1 

K-S Test Results'^ 



o 

II 

Tb = 0.0 — 0.3 

Tb = 0.3 - 1.0 

/^phot 

0.06 

0.008 

0.09 

dSED 

0.05 

0.004 

0.05 


^ Probability that the distributions of /3 in 
the two bins of log(sSFR/yr“^) (—9.60 < 

log(sSFR/yr“^) < —8.84 and —8.84 < 

log(sSFR/yr“^) < —8.00) in the given bins of Tb 
are drawn from the same parent population. 

older and more massive. For the moment, we consider 
the effect of star-formation history (i.e., effect of older 
stellar populations) on /3 by color-coding the points in 
Figure by sSFR. It is apparent that there is a “se- 
quence’^in the location of galaxies in the Pphot-Th plane 
with sSFR, such that the locus of points shifts towards 
redder ,5phot at a fixed Tb with decreasing sSFR. This is 
the expected behavior if late-type stars increasingly af¬ 
fect the near-UV flux with decreasing sSFR. The shifting 
of the locus is also (unsurprisingly) apparent if we differ¬ 
entiate the points again by sSFR, but use /3sed instead 
of /Jphot (Figure!^. 

To demonstrate quantitatively the shifting of the /? vs. 
Tb relation with sSFR, we executed two tests. In the first, 
we used K-S statistical tests to determine the probability 
that two distributions of /3 at two given sSFRs (and at a 
fixed Tb) are drawn from the same parent distribution. In 
practice, we gridded the data in three bins of Tb (Tb=0, 
0.0-0.3, and 0.3-1.0), and two bins of log(sSFR/yr“^), 
where —9.60 < log(sSFR/yr“^) < —8.84 and —8.84 < 
log(sSFR/yr“^) < —8.00. The sSFR bin boundaries 
were chosen to encompass most of the objects in our 
sample, with an equal number of objects per sSFR bin. 
We then performed two-sided K-S tests on the two dis¬ 
tributions of P (for the two sSFR bin^ in each bin of Tb. 
The results are summarized in Table and they suggest 
that the P distributions between the two sSFR bins are 
drawn from different parent populations. 

In the second test, we divided the data into the same 
two bins of sSFR defined above, and we computed 
the Spearman correlation coefficient and the significance 
with which the null hypothesis of no correlation between 
P and Tb can be ruled out. Table summarizes the re¬ 
sults for /3sed vs. Tb and /3phot vs. Tb. Significant correla¬ 
tions between P and Tb are found in most cases, with the 
slopes and intercepts from linear regressions indicated in 
Tablel^ We also estimated the intrinsic dispersion in the 
data (CTint) by assuming that the total dispersion in P vs. 
Tb is the quadrature sum of dint and the measurement 
error, Umeas- Clearly, dint will depend on the adopted 
bin widths, but our goal was to determine the degree to 
which the P — Th correlations in given sSFR bins overlap 
when considering the intrinsic scatter in the data. The 
fits obtained for /3sed vs. Tb are shown in Figurej^ From 
these statistical tests, it is clear that the scatter seen in 
Figure is to a large extent due to variation in the P 
vs. Tb rdations for galaxies with different sSFRs. Specif¬ 
ically, the relation for the lowest sSFR galaxies is shifted 
to redder /3 by « 2 x dint relative to the relation for the 
highest sSFR galaxies. 

3.3.2. Contribution of A Stars to the Near-UV Flux 
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Figure 2. Photometric {left-, /3phot) SED-based {right, /3sed) UV slopes vs. Balmer optical depth (ri,), differentiated by sSFR. Errors 
in Tb include both measurement uncertainties in the Ha and H/3 fluxes, and the uncertainty in slitloss corrections. Errors in ftphot are 
derived directly from the linear regression between log fx and A. Errors in /3sed are estimated from Monte Carlo simulations of the SED 
fitting (see text). The sSFRs assume dust-corrected Ha-based SFRs and the stellar masses returned from the SED fitting. The large stars 
indicate the average positions of the H/3-undetected galaxies. Galaxies with rn = 0 and different ranges of sSFR are arbitrarily offset to 
the left side of each panel for clarity. 


Table 2 

Spearman Correlation Test and Linear Regression Results for /3 vs. ru 


0 

log[sSFR/yr-l] 

Ara 

PSpear 

^Spear^ 

Slope 

Intercept 



/^phot 

-9.60...-8.84 

87 

0.33 

3.12 

0.57 ± 0.15 

-1.42 ± 0.02 

0.52 

0.33 

-8.84...-8.00 

87 

0.56 

5.24 

0.73 ±0.11 

-1.67 ± 0.03 

0.38 

0.21 

dsED 

-9.60...-8.84 

87 

0.55 

4.42 

0.95 ±0.14 

-1.48 ±0.02 

0.47 

0.31 


-8.84...-8.00 

87 

0.65 

6.07 

0.87 ± 0.09 

-1.78 ±0.03 

0.33 

0.20 


^ Number of galaxies. 

Spearman rank correlation coefficient. 
Standard deviations from null hypothesis, 
d Total RMS. 

® Intrinsic dispersion (see text). 


We have shown that galaxies with lower sSFRs have 
on average redder UV slopes at a given Balmer optical 
depth, but it remains to be seen whether late-type stars 
are responsible for this trend. To assess the contribution 
of lower mass A stars to the UV continuum, we recalcu¬ 
lated the UV slopes using only those photometric points 
that lie blueward of rest-frame 1750 A, /3(1260— 1750 A). 
Only O and B stars, and the most massive A stars (spec¬ 
tral type Al, AO) contribute significantly at these blue 
wavelengths (e.g., iLeitherer & Heckman 119951 ICalzetti 


et al. 1994). We compared these “blue' UV slopes to 


those obtained using photometry over the full wavelength 


range between 1260 and 2600 A, /3(1260 — 2600 A). The 
comparison (Figure 13 includes only those galaxies that 
had a minimum of 4 photometric points available for cal¬ 
culating /? when excluding all points with A > 1750 An 
We find no evidence that galaxies with lower ss- 


FRs have /lphot(1260 — 1750 A) bluer than /?phot(1260 — 
2600 A), as would be expected if the near-UV continua 
of such galaxies were weighted more heavily by the flux 


® The galaxies shown in Figured are all in the COSMOS field, 
where 0 is determined using botnbroad- and intermediate-band 
photometry with sufficient independent wavelength sampling to ac¬ 
curately compute the UV slope in the relevant wavelength ranges. 


from A stars. The same conclusion is reached when con¬ 
sidering the relation between /?phot(1260 — 1750 A) and 
the slope obtained using photometry strictly in the near 
UV, /lphot(1750 - 2600A). 

The redness of the UV slope measured over the shorter 
baseline in wavelength can in principle reflect the more 
severe attenuation at these short wav elengths, an issue 
that we discuss briefly in Section |4.2| Because the dif¬ 
ference between (3 measured across the short and long 
wavelength baselines does not appear to correlate with 
sSFR (Figure]^, the contribution of older stars to the 
near-UV continuum is likely minimal over the range of 
sSFR considered here, and thus not a dominant factor in 
the spread or trend in /3 with rb. 


3.4. Effect of Metallicity on the UV Slope 

An issue separate from the effect of older stars on the 
near-UV continuum is the extent to which differences in 
metallicities of the stellar populations lead to trends be¬ 
tween P and Tb. We investigated the dependence of /3 
on metallicity by using MOSDEF measurem ents of the 
[NII]A6585/H a ratio (“N2 index”) with the Pettini ^ 


Pagel (2004) calibration to infer the gas-phase oxygen 
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Figure 3. Linear regressions between /3 sed The solid 

lines and shaded regions indicate the best-fit linear functions and 
estimated zhcint intrinsic scatter for the different bins of sSFR, 
respectively. Individual data points have been omitted for clarity. 
Long-dashed lines indicat e the best-fit l i near f unction and upper 
and lower envelopes from |Calzetti et al.| (|1 994|| for a sample o f 39 
local starburst and blue compact galaxies i jKinney et al.|1993[ |. 



/3p,„t(1260-2600A) 


Figure 4. Comparison of UV spectral slopes computed over dif¬ 
ferent wavelength ranges that include (/3phot(1260 — 2600 A)) and 
exclude (/3phot(1260 — 1750A)) the near-UV portion of the spec¬ 
trum that can have a significant contribution from stars of A2 and 
later spectral type. We find no evidence that /3phot(1260 — 1750 A) 
is systematically bluer than ^phot(1260 — 2600 A), suggesting that 
A stars do not contribute significantly to the near-UV continuum, 
even for galaxies with lower sSFRs. 


abundance 13 We have assumed that the gas-phase and 
stellar met^icities track each other closely, as should be 
the case for the short-lived O and B stars that dom¬ 
inate the UV continuum. The range of metallicities 
probed by our sample is 8.1 <12-1- \og{0/H) < 8.7 
(i.e., 0.2 < Z < 1. 0ZfT>, assuming the value o f the so- 
lar abundance from Allende Prieto et al. (2001|) and As- 


plund et al. (2004)), with a mean oxygen abundance that 


changes by ~ 0.3 dex between the reddest and bluest 


® A discussion of the differen t metallicitv indicato rs accessible 
with MOSDEF are discussed in|Sanders et al.| (|2015|l. Assuming 
these different indicators does not alter our conclusions regarding 
the change in 0 with inferred oxygen abundance. 


UV slopes. For a 

Chabrier 

(2003) IMF and a constant 

star formation for 

lUU Myr, the 

Bruzual & Chariot 1 

2003) 


stellar population synthesis models imply a ditterence in 
intrinsic UV slope of A/3int — 0.2 for the range of metal- 
licites probed in our sample. This difference is much 
smaller than the spread in /3 at a given Tb, and much 
smaller than the change in /3 as Tb increases (Figure [3). 
We conclude that differences in the metallicities of the 
stellar populations are unlikely to be the dominant fac¬ 
tor in driving the change in UV slope with increasing 
Balmer line opacity. 


3.5. Scatter in /? in the limit where Tb —>■ 0 

The previous discussion of the effects of late-type stars 
and metallicity on the UV slope leads to the conclusion 
that dust obscuration is likely the underlying cause of 
the trend between /3 and Balmer optical depth. In this 
scheme, variations in the dust opacity between individ¬ 
ual galaxies are likely responsible for the scatter in /3 at 
a fixed Tb. Such variations in the opacity may be due 
to a dependence on dust composition, resulting in differ¬ 
ent absorption and scattering cross-sections of the dust 
grains. Alternatively, the geometrical distribution of the 
dust with respect to the ionizing stars in galaxies can 
lead to variations in dust opacity. Galaxies of a given Tb 
may exhibit a wide range of /3 depending on the dust col¬ 
umn densities to wards the ionized regions a nd the stellar 


continuum (e.g., Calzetti et al. 


1994 2000). Though, as 
emphasized in |(Jalzetti et al. Il994p, the spatial distri¬ 
bution of dust and stars cannot explain the large scatter 


explain the large i 
in j3 in the limit where Tb — >■ 0. In this limit, we ex¬ 
pect /3 to approach its intrinsic value, typically between 
/3int — —2.6 and —2.0, depending on the details (e.g., 
IMF, metallicity) of the stellar population model. 

There are several reasons to suspect that galaxies with 
Tb = 0 are not dust-free, even though their Balmer decre¬ 
ments would indicate formally otherwise. First, such 
galaxies have measurement errors in Tb that do not ex¬ 
clude the presence of some dust. Second, even galaxies 
with the bluest UV spectral indices in our sample have 
metallicities of at least 12 -\-\og{0/H) > 8.1 (« 0.26 Zq) 
implying that they must contain some dust. Third, if we 
de-red den the best-fit SEDs assuming a [Calzetti et al.| 
(20001 attenuation curve and the best-fit E{B — V), we 
obtain /3int — —2.3. This is significantly bluer than the 
Ant ~ —1.3 and —1.6 implied by the linear regression for 
the bin of lowest sSFR (Tablej^). Moreover, the detection 
of Ha emission in these galaxies requires the presence of 
hydrogen-ionizing stars, and the blue hot stars will cause 
an otherwise intrinsically shallow UV spectral index to 
become steep (i.e., more negative). Finally, from an ob¬ 
servational standpoint, if the rest-UV emission in these 
galaxies is not spatially coincident with the rest-optical 
emission, then the observed Balmer decrement may not 
be indicative of the global reddening of the galaxy as 
deduced from the UV slope, a case that may apply for 
all galaxies in our sample. For these reasons, the linear 
regressions between /3 and Tb should be interpreted with 
caution in the limit Tb —>■ 0. 


To_sumimrize, the results of Section 3.3.1 (e.g., Fig¬ 
ures an d§ imply that the physical origin of the scatter 
in P vs. Tb must vary systematically with sSFR, a point 
that we discuss further in the next section. We have 
investigated the potential effect of A stars on our mea- 
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surements of /3 and find that it cannot explain the scatter 
in /3 at a given Further, variations of the high-mass 
end of the IMF (e.g., as a function of m etallicity) lead 
to co rrections of the UV slope of « 10% (Calzetti et al. 
19941, which again cannot account for the large scatter 
in d at a given rt- In Section]^ we put forth a physi¬ 
cal interpretation that is consistent with the trends and 
scatter between /3 and Tb- 

3.6. The Attenuation Curve of High-Redshift Galaxies 

Given evidence that neither an older stellar popula¬ 
tion nor a change in metallicity is responsible for the 
reddening of (3 with increasing Tb, we are left with the 
possibility that dust attenuation is the dominant factor 
in the observed trend be tween and Tb . A sim ilar line of 
reasoning is presented in Calzetti et al. (1994) to account 
for the correlation between U V slope and Daimer optical 
depth for a sa mple of 39 local starb urst and blue compact 
galaxies from Kinney et al. ( 1993|), where the local trend 
is shown in Figure!^ Interestingly, we find that the cor¬ 
relation between UV slope and Balmer optical depth is 
steeper for the local sample, relative to that observed for 
the 1.4 ^ z < 2.6 MOSDFF sample. The median sSFR of 
the local galaxies is log(sSFR/yr“^) « —9.00, and these 
galaxies generally have both lower SFRs (~ 4M0yr“^) 
and lower stellar masses (~ 3.9 x 10® Mq) compared to 
those in our sample. A more detailed comparison of our 
results with local samples (including the SDSS) will be 
presented elsewhere. For the moment, we note that the 
segregation of the /3 vs. Tb relations (Figure]^ indicates 
an sSFR-dependent scaling between nebular and stellar 
attenuation in hi gh-red shift galaxies that differs from lo¬ 
cal ones (Section 3.6.2). 


3.6.1. Derivation of the Selective Attenuation Curve 

In what f ollows, we adopt a pr ocedure similar to that 
p resented in Calzetti et al.| (19941, generalized from that 
of Kinney et al. ] 1994), to derive the selective attenuation 
curve. 'I'he premise behind this procedure is that because 
the UV spectral index is sensitive primarily to dust at¬ 
tenuation, the ratio of the average spectrum of galaxies 
with high Balmer optical depths to that of galaxies with 
low optical depths should reflect the average effect of at¬ 
tenuation on the stellar SEDs. Of course, especially in 
the optical and near-IR, there can be a wide range of 
intrinsic SED shapes for star-forming galaxies depend¬ 
ing on their star-formation histories and ages. How¬ 
ever, to first order, galaxies of a given sSFR will have 
comparably-shaped intrinsic stellar SEDs, as the flux 
at any given wavelength is contributed to in a similar 
proportion by the current star formation and the light 
from previous generations of star formation. The ex¬ 
act form of how the previous stellar mass was built up 
(i.e., the star-formation history) is unimportant, as the 
rest-frame optical to near-IR light that is used to infer 
stellar masses (and hence sSFRs) arises from some com¬ 
bination of the currently forming stellar population, and 
from lower mass (spectral type F5 and later) stars that 
have not yet evolved off the main sequence for the age of 
the Universe at z ~ 2{^ For instance. Figure [s] shows the 

This assumes an IMF that does not evolve with redshift above 
2 ~ 2. Some circumstantial evidence for an un-evolving IMF at 



Figure 5. Stellar SEDs for a galaxy with a stellar mass of 
M* ~ 1.2 X 10^'^ M 0 and SFR~ 16 MQyr~^, assuming differ¬ 
ent simple star-formation histories; an exponentially rising model 
with an e-folding time of 5000 Myr, constant star formation, and 
exponentially declining models with e-folding times of 1000 and 
100 Myr. The corresponding ages of the stellar populations are Ri 
800, 3000, 720, 570, and 200 Myr, respectively. The SEDs are cast 
in terms of fi, and are normalized so that A (A = 5500 A) = 1 . 
Fixing the sSFR of the galaxy results in similar SED shapes for 
most of the assumed star-formation histories. 

conformity of the stellar SED for a galaxy with an SFR 
and M* (and hence, sSFR) typical of the average galaxy 
in our sample, irrespective of the star-formation history. 
SEDs that assume a relatively short decay timescale for 
star formation (e.g., Td < 100Myr) deviate by « 10% at 
optical and near-IR wavelengths, but such histories are 
ruled out for the vast majority of objects in our sample 
given their typical stellar population ages of > 300 Myr, 
combined with the evidence that, at least on average, 
> L* galaxies at z ~ 2 show evi dence for rising sta r 
formation at previous epochs (e.g., Reddy et al.||2012b). 

Motivated by the apparent sSFR-dependence ot the 
j3 vs. Tb relation, we derived the attenuation curve for 
galaxies in two bins of sSFR. We first examined the dis¬ 
tribution of sSFR with Tb (Figure!^. Approximately 
96% of the galaxies in our sample (174 of 182) lie in 
the range —9.6 < log[sSFR/yr“^] < —8.0. Dividing the 
number of objects in this range into two bins, with equal 
numbers of objects (87) per bin, results in the following 
sSFR bins: -9.60 < log[sSFR/yr-i] < -8.84 (sSFR bin 
1) and -8.84 < log[sSFR/yr-i] < -8.00 (sSFR bin 2). 
The objects within each sSFR bin were further divided 
into three bins of Tb. For sSFR bin 1, the first Tb bin in¬ 
cludes all objects with Tb = 0, and the remaining objects 
are split roughly equally into two additional bins. Given 
the low number of objects (8) that have Tb = 0 in sSFR 
bin 2, for this sSFR bin we divided the objects equally 
into three bins of Tb. 

The derivation of the attenuation curve relies on the 
assumption that the average intrinsic SED (and hence 
the sSFR) is unchanging as a function of Balmer opti¬ 
cal depth, irrespective of the exact mix of different star- 
formation histories that may be contributing to that av¬ 
erage SED. Figure shows that the average sSFRs of 
galaxies in each bin of Tb is roughly constant, imply¬ 
ing that a simple ratio of the average SEDs obtained in 

2 > 2 comes from the general agreement between the integral of the 
cosmic SFR density and the evolution of the stell ar mass density 
( [Reddy & Steidel|2009l |Madau fc Dickinson|2014[ |. 
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Figure 6. Distribution of sSFR with Balmer optical depth. Ob¬ 
jects with Tb = 0 are offset to the left for clarity. The horizontal 
lines denote the boundaries of the two sSFR bins and the vertical 
lines denote the boundaries of the three rt, bins per sSFR bin used 
to calculate the attenuation curve. The open diamonds indicate 
the mean values of log[sSFR/yr~^] and Tb for each bin. 


each bin will directly probe the wavelength dependence of 
dust obscuration. In Section |3. 6 . 4[ we address the slight 
(< 0.1 dex) differences in the (log[sSFR/yr“^]) obtained 
in each bin, and how these differences translate into a 
systematic uncertainty in our derived attenuation curve. 

Having established the aforementioned bin definitions, 
we then estimated the average SED of objects in each 
bin as follows. The broad- and intermediate-band fluxes 
measured for each object were normalized by the flux of 
the best-fit SED for that object at rest-frame A = 5500 A 
and then shifted to the rest frame based on its spectro¬ 
scopic redshift. The normalized flux density measure- 
ments for all object s were then combined and fit with a 


Bruzual & Chariot (2003) stellar population model. In 
the htting, we assumed an equal weighting for all ob¬ 
jects irrespective of the flux measurement errors so that 
no single object dominated the fit, though we note that 
taking into account photometric errors results in fits that 
are similar to those obtained here. To ensure the best 
possible fit to the data, we allowed for a range of pos¬ 
sible star-formation histories (declining, constant, and 
rising), no restriction on the age of the model apart from 
it being less than the age of the Universe at the mean 
redshift of the sample {z) ~ 2 . 0 , color excesses from 
E{B — V) = 0.00 to 0.60, and four attenuation curves, a 
Calzetti starburst reddening curve (Calzetti et al.|2000), 
and the line-of-sight SMC, EMC, and Galactic extinction 
curves. The exact details of the best-fit model are not 
essential, so long as the model provides a good fit to the 
fluxes of the contributing objects. Figure shows the 
best-fit functional form to the normalized photometry 
of all objects in each bin of sSFR and Tb. In this re¬ 
gard, also shown are the weighted mean fluxes of objects 
in 1000 A- width bins in wavelength (i.e., the empirical 
SED), demonstrating that the model fits are not in any 
way skewed relative to these mean empirical values. Be¬ 
low, we also demonstrate that our model fits result in 
a shape of the attenuation curve that follows the one 
calculated from the empirical SEDs. 

The random uncertainty in the average SEDs was es¬ 
timated as follows. For each object, we perturbed its 
sSFR and Tb based on the errors in these quantities. 


and then rebinned all the objects according to these per¬ 
turbed values. The photometry of each object was also 
perturbed according to the photometric errors, and the 
average SEDs were recalculated. This procedure was re¬ 
peated 100 times to obtain statistically robust determi¬ 
nations of the standard deviations in the average SEDs, 
which are shown by the shaded regions in Figure 
The selective attenuation between each of these average 
(composite) SEDs is then defined as 


Q^{X) = — 


Ti - Ti 


In 


m) 


FiiX)\ 


( 2 ) 


where Ti and Ej refer to the average Balmer optical depth 
and the composite spectrum, respective ly, of the objects 
in the ith bin of Tb (Calzetti et al.|I994|. Qi is computed 
relative to the first bin where by construction ti = 0 for 
the first sSFR bin, and ti = 0.065 for the second sSFR 
bin. Thus, for each bin of sSFR, we obtain two esti¬ 
mates of the selective attenuation (Q 2 and Q 3 ). Figure 
shows our calculation of Q 2 and Qa for the two bins a 
sSFR specified above. The uncertainties in Qi were de¬ 
termined by randomly perturbing the composite SEDs 
according to their respective uncertainties and recom¬ 
puting Qi for these realizations. The errors in the mean 
Qi are indicated by the shaded regions in Figure and 
are cr{Q 2 ) « 0.20 and ^(Qa) « 0.05 for sSFR bin I, and 
o'{Q 2 ) ~ O.IO and cr^Q^) « 0.04 for sSFR bin 2. The 
differences between Q2 and Qa for a given sSFR bin are 
^Qi ^ 0.15, or < 0.07 dex, and are of the same order as 
the random uncertainties in the mean Qi. 

Additionally, we calculated the Qi based on the em¬ 
pirically determined SEDs. This was accomplished by 
linearly interpolating between the average flux points of 
each composite SED (i.e., as shown by the diamonds in 
Figure]^ to produce an empirical SED. These empirical 
SEDs were then combined in the same way as the model 
fits in order to calculate Qi (i.e., using Equation!^, and 
the results are denoted by the dashed lines in the left 
panel of Figure (for clarity, we show the results for 
Qs). The empirically determined Qi track closely those 
inferred from our model fits, and for our subsequent anal¬ 
ysis, we adopt the latter. 

Our analysis indicates that the Qi within a given sSFR 
bin are very similar in shape, implying that the average 
SEDs used to compute Qi have roughly similar stellar 
populations (as would be expected given the discussion 
above regarding the similarity in SED shapes for a given 
sSFR). In other words, the shape of Qi is dictated pri¬ 
marily by differences in the dust obscuration, as param¬ 
eterized by Tb, between the average SEDs show in Fig¬ 
ure and not by significant differences in the average 
intrinsic stel lar po pulations from one Tb bin to another 
(see Section 3.6.4 for a full discussion of the systematic 
uncertainties in tne attenuation curve). Given the sim¬ 
ilarity in the shapes of Q 2 and Q 3 for each sSFR bin, 
we simply took a weighted average of the two selective 
attenuation curves to produce what we refer to as the 
“effective” attenuation curve, QeS- The weights were cal¬ 
culated as I/(t|, where (T\ is the uncertainty in Qi as a 
function of wavelength. The uncertainty in the weighted 


average is cTweighted.A = ^ 1 / Z) o'a 


We fit the effective attenuation curve in each sSFR 






















MOSDEF: The Dust Attenuation Curve at z ~ 2 


11 



Figure 7. Photometry (circles) and model fits (red) of objects in bins of sSFR and Balmer optical depth, with number of contributing 
objects and mean redshifts indicated. The photometry has been normalized such that /;^ (5500 A) = 1 . The fitting takes into account the 
photometric (flux) errors and allows for a range of possible star-formation histories, ages, color excesses, and attenuation curves (see text). 
The orange shaded regions indicate the 1 cr dispersion about the mean SED fits, based on Monte Carlo simulations where the sSFRs, rb, and 
photometric errors for each object are perturbed and the mean SED is re-calculated. For comparison, the diamonds show the un-weighted 
average of the fluxes in 1000 A-width bins in wavelength. 




Figure 8 . Left: Selective attenuation curves {Q 2 and Q 3 ) calculated using objects in two bins of sSFR. The curves for sSFR bin 2 have 
been shifted upward by = 1 for clarity. The shaded regions indicate the uncertainties in Qi. The Qi inferred from the empirical SEDs 
are denoted by the dashed lines. Also shown at the bottom is the difference between Q 3 and Q 2 as a function of wavelength for each bin 
of sSFR. Right: Effective attenuation curves (Qeff) calculated as the weighted average of the Q 2 and Q 3 s hown in the lef t pane l for each 
sSFR bin. The weighted uncertainties in QefF ^re indicated by the shaded regions. The dashed line shows the|Calzetti et al.|(|1994|| selective 
attenuation curve. 


bin using third-order polynomial functions of 1/A over 
the wavelength ranges A = 0.15 — 0.60//m and A = 
0.60 — 2.20 ^m. These specific ranges were chosen to 
ensure the best possible representation of the data using 
polynomials of the smallest order. The polynomial fit for 
sSFR bin 1 (-9.60 < log[sSFR/yr-i] < -8.84) is: 

Qeff(A,-9.60 < log[sSFR/yr-i] < -8.84) 

= -2.101 + 1.450/A - 0.181/A2 + 0.010/A^ 

0.15 < A < 0.60 /.tm; 

= -0.955 - 0.169/A + 0.732/A^ - 0.200/A^ 

0.60 < A < 2.20 Aim. (3) 

The polynomial fit for sSFR bin 2 (—8.84 < 


log[sSFR/yr-i] < -8.00) is: 

Qeff(A,-8.84 < log[sSFR/yr-i] < -8.00) 

= -1.834 + 1.299/A - 0.179/A^ + 0.010/A^ 

0.15 < A < 0.60 Atm; 

= -0.878 + 0.136/A + O.M8/X'^ - 0.091/A^ 

0.60 < A < 2.20 Atm. (4) 

The error in Qeff at each wavelength point was used to 
weight the data appropriately in the fitting. The formal 
uncertainty in the model fit at each wavelength varies 
in the range ctfit « 0.001 to 0.004. As discussed be¬ 
low, several additional steps are required to convert QeS 
to a total attenuation curve, the reason being that the 
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former is computed using galaxies of different luminosi¬ 
ties at different redshifts, and therefore the normalization 
of the attenuation curve must be determined separately. 
Nonetheless, Qeff is useful in two respects: it indicates 
the wavelength dependence of the obscuration of the stel¬ 
lar continuum and the optical depth to dust of the con¬ 
tinuum relative to the Balmer lines. We discuss these 
points below. 


3.6.2. Comparison of the Average Optical Depths of the 
Continuum and the Balmer Lines 

The definition of Q^s (Equation implies that the 
Balmer optical depth, Tb, is related To the difference in 
optical depths of the continuum at the wavelengths of H/3 
and Ha: rr.(Hd)-Tr( HQ) = [Qeff (H/3) - Qeff (Ha)]Tb (e.g., 
Calzetti et al.||1994), assuming that (rb) = 0 for the first 
bin of Tb. Accordiiig to Equations and we can then 
write Tc(H/ 3) — T(.(Ha) = (0.42 ± 0.01)rb and « (0.36 ± 
0.01)Tb for sSFR bins -9.60 < log[sSFR/yr-i] < -8.84 
and —8.84 < log[sSFR/yr“^] < —8.00, respectively]^ 
Thus, we conclude that analogous to loca l galaxies (e.g., 
Fanelli et al .j|T9 88j [Calzetti et al. 1994[ |Mas-Hesse fc| 
Kreckel et al.|2013D, the ionized gas in high- 


Kunth|1999[ 


. gas in high- 

redshift galaxies with stellar masses > 10® Mq is subject 
to a higher dust optical depth on average than the stellar 
continuum. Qualitatively, this result is consistent with 
the average difference in color excess between the gas 


and stars observed in local star- forming galaxies (Calzetti 
et al.|[2000 Kreckel et al.||20i3 ). We return to this point 


in Section |4.3[ where we investigate in detail the rela¬ 
tionship between the total attenuation of the ionized gas 
and stellar continuum for individual galaxies. 


3.6.3. Shapes of the Attenuation Curves 
The total attenuation curve is defined as 


fc(A) = 


Aa 

E{B - V) ’ 


(5) 


where A\ is the total magnitude of extinction at wave¬ 
length A and E{B — V) is the color excess of the stellar 
continuum. Because E{B — V) = A{B) — A(V), we can 
write k{B) — k{V) = 1. Thus, we multiplied the effec¬ 
tive attenuation curves Qeff by the factors / = 2.676 and 
3.178 for sSFR bins 1 and 2, respectively, so that 


f[QMB) - geff(E)] = k{B) - k{V) = 1. (6) 


The factors were computed by evaluating the effective 
attenuation curves in the B {X = 4400 A) and E (A = 
5500 A) bands. The resulting /Qeff are shown for the 
two bins of sSFR in Figure]^ 

The two notable aspects of the effective attenuation 
curves derived here are that (a) they a ppear to devi¬ 
ate at long wavelengths relative to the [Calzetti et al.| 
(|1994|) selective attenuation curve, and (b) QeS for the 
two sSFR bins are strikingly similar. On the latter point, 
a priori we would not have expected such a good agree¬ 
ment given the disparity in the trend of UV slope with 
Balmer optica l depth for galaxies with different sS FRs 


(Section 3.3.1). We return to this issue in Section 4.3 


The relation is only approximate for bin where —8.84 < 
log[sSFR/yr“^] < —8.00 as the average Balmer optical depth for 
the bin of lowest rb is not exactly zero. 



Figure 9. fQeS for tho two bins of sSFR as indicated in the 
figure, where the multiplicative factors / = 2.676 and 3.178 are 
required to ens ure that f\Q^fr(B) — Q ^fr(V)] = k{B) — k(V) = 1 
(see text). The [Calzetti et al.|||1994|| attenuation curve is shown 
for comparison, whe re we nave applied a similar factor / Ri 2.659 
i jCalzetti et al.|2000| |. 

Given the similarity in Qeff derived for the two sSFR 
bins, we simply averaged them, resulting in the following 
functional form for the total attenuation curve: 


k{X) = -5.726 -b 4.004/A - 0.525/A® -b 0.029/A® -b Rv, 
0.15 < A < 0.60/xm; 

= -2.672 - 0.010/A -b 1.532/A® - 0.412/A® -b Rv, 
0.60 < A < 2.20/xm; (7) 


where Ry has the usual definition as the ratio of total- 
to-selective absorption at E-band: Ry = Ay/E{B — V). 
The last step in the derivation of the total attenuation 
curve is to calculate Ry . This can be accomplished by ex¬ 
trapolating k{X) to some wavelength (e.g. Xn = 2.85 /xm) 
and assuming that k(X > Xn) ~ 0. This particular value 
of Xn = 2.85 /xm corresponds to where k(X) for t he SMC, 
EMC, and Milky Way are very close to zero ( jCor donj 
et al.|200'3). In our case, if we extrapolate EquatrnnjTjto 


Xn = 2.85 /im and force fc(A„) = 0, then Ry = 2.505~ 

Alternatively, a common second approach is to extrap¬ 
olate k(X) to infinite wavelength, and set k(X —>■ oo) = 0. 
Doing so for Equation results in Ry = 2.672. More 
typically, a k(X >> 1 /xmj ex A“® dependence is generally 
assumed. In this case, we determined the functional form 
of k(X > 2.20 ^m) by fitting Equation^with a power-law 
in 1/A in the range A = 1.9—2.2 /xm (this is to ensure that 
k(X) is a continuous function across A = 2.2/xm). With 
this long-wavelength extrapolation of fc(A), and setting 
k(X —>■ oo) = 0, we find Ry = 2.942. 

The different extrapolations discussed abov e re sult in 
a systematic error in Ry of « 0.4. Figure shows 
the comparison of the attenuation curves assuming the 
smallest and largest values of Ry derived above, along 
with several other common extinction and attenuation 
curves. Independent estimates of Ry based on mid- and 
far-infrared data will be presented elsewhere. For the 
remaining analysis, we have assumed the first value of 
Ry = 2.505 as it does not require an extrapolation too 
far beyond where we have constraints on the shape of the 
attenuation curve from the photometry (A ^ 2.2/xm). 
With this assumption, we arrive at the final expression 
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Figure 10. Left: Compar i son of the MOSDEF attenuation curve assuming two different values of Ry = 2.505 and Ry = 2.942 (see text), 
and that of|Calzetti et al.| ||2000|l {dotted line). Also inclu ded are the local ext inction curves of the Milky Way, SMC, and LMC, where 
the shapes and normalizations [Ry) have been taken from|Gordon et al.| (|2003||. Right: Same as left panel, where we show the MOSDEF 
attenuation curve assuming Ry = 2.505 over the full wavelengtn range A = U.ii — 2.50 ^m. 


for the total attenuation curve relevant for the stellar 
continuum: 


fcstars(A) = -5.726 + 4.004/A - 
+0.029/A3 + 2.505, 

0.15 < A < 0.60/xm; 

= -2.672 - 0.010/A + 1.532/A2 
-0.412/A3 + 2.505, 

0.60 < A < 2.85 ^m. (8) 

3 . 6 . 4 . Systematic Uncertainty in k{X) Arising from 
Differences in Average Stellar Populations 


As noted previously, the procedure we have used to de¬ 
rive the attenuation curve assumes that the average in¬ 
trinsic SEDs in the different bins of Balmer optical depth 
do not vary substantially from one another. Our method 
of binning the galaxies by sSFR was designed to mitigate 
any such intrinsic spectral variations. However, while the 
mean sSFRs as a f unction of rb are similar from bin-to- 
bin (Section 3.6.1 and Figure [^, they are not identical 
and will therefore introduce some bias in our estimate of 
the attenuation curve. 

To quantify this bias, we first computed an “intrin¬ 
sic” SED for each galaxy by: (a) assuming a rising star- 
formation history with an e-folding time of Tr = 5000 Myr 
(this is essentially equivalent to assuming a constant 
star-formation history); (b) fixing the stellar mass of 
the stellar populatio n to be the value derived from the 
SED-fitting (Section [2^; and (c) varying the age of the 
stellar population until the corresponding SFR of the 
model is equal to the (B alm er decrement corrected) Ha- 
determined one (Section[2^. The goal of these steps was 
to construct an SED for a stellar population that has the 
same sSFR as that computed from SFR(Ha) and M*. 
Hence, the exact form of the star-formation history as¬ 
sumed in step (a) above is inconsequential, as per the 
discussion in Section r3. 6 . II 

Second, these intrinsic SEDs were normalized so that 
/(A = 5500 A) = 1 and then averaged together for all 
objects within each bin of sSFR and Tb (assuming the 
sa me b ins used to derive the attenuation curve). Fig¬ 
ure [m shows these average intrinsic SEDs. The differ¬ 
ences in these intrinsic SEDs lead to “correction factors” 


AQint that must be applied to the Qi derived above: 
AQint = l./(r, - Ti) X ln(/”V/“), where is the 
intrinsic spectrum inferred for the A** bin of Tb. These 
corrections factors, as a function of wavelength, are also 
show in the bottom of the left panel of Figure El 

Third, we computed the corrected Qi as Q™*' = Qi + 
AQint- We then followed the same procedure discussed 
above to compute fQlg- i.e., the Q^™ in each bin of 
sSFR were averaged together to produce for each 
sSFR bin, third-order polynomials were fit to <3“^ , and 
these functions were multiplied by the relevant correc¬ 
tion factors according to EquationFigure [^indicates 
that the fQl^ are essentially identical to the uncorrected 
versions initially shown in Figure In effect, correcting 
for the small differences in the intrinsic SEDs is mostly 
compensated for by the multiplicative factors / which 
are smaller for the corrected Qi. The ratios of the un¬ 
corrected to corrected fQeS at A = 0.15 and 2.20/xm are 
1.083 and 0.998, respectively, for sSFR bin 1; and 0.992 
and 1.056, respectively, for sSFR bin 2. Given the less 
than « 8% differences implied for the total attenuation 
curve, we assume the uncorrected version expressed in 
Equation 

Finally, we note that systematic offsets in the derived 
stellar masses or SFR(Ha) resulting from the assump¬ 
tion of different attenuation curves in the SED fitting 
will simply shift all the galaxies up or down in terms of 
sSFR, but the ordering of the galaxies will for the most 
part be preserved!^ For this reason, the ratios of the 
intrinsic SEDs shown in Figure E] will be unaffected by 
assuming a different attenuation curve in the derivation 
of SFR(Ha) and M*. Detailed comparisons of the SED 
parameters and SFRs derived with different assumptions 
of the attenuation curve are presented in Section 


3 . 6 . 5 . Comparison to Other Common Curves 

The average attenuation curve found for the MOSDEF 
galaxies is compared to several local l ine-of -sight extinc¬ 
tion curves and the Calzetti et al. (2000) attenuation 


The ordering of galaxies will be preserved if the same atten¬ 
uation curve applies to all galaxies. This is the case for galaxies 
in our sample given the similarity in shape between Q 2 and Q 3 in 
a given sSFR bin (Figure]^, and the similarly in shape of /Qeff 
between the two sSFR bins (Figure |^. 
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Figure 11. Left: Inferred average intrinsic SEDs for the same bins of sSFR and rb used to calculate the attenuation curve. The intrinsic 
SED for each galaxy is constructed by assuming a simple stellar population with the same SFR(Hq) and M* as measured for that galaxy. 
These individual SEDs of objects in each bin of sSFR and rt, are then averaged together to produce the curves shown in the Figure. The 
bottom panel shows the implied correction factors to the selective attenuation curves due to differences in the intrinsic SEDs of galaxies 
from bin-to-bin. Right: Comparison of the /Qeff {thin lines) with the versions corrected for variations in the intrinsic SEDs {thick lines), 
for each bin of sSFR. 


curve in Figure The SMC, LMC, and Milky Way 
curv es, and their respect ive normalizations, are taken 
from Gordon et al. (2003). A common feature of these 
curves is that they all asymptote to zero as the dust ab¬ 
sorption cross-section becomes negligible at long wave¬ 
lengths, with a corresponding rapid rise in the attenu¬ 
ation curve toward shorter (UV) wavelengths. The one 
distinguishing feature of the Milky Way and LMC extinc¬ 
tion curves, namely the 2175 A “b ump” (Stecher & Bonn 
1965| |Fitzpatrick & Massa||1986|) , is absent from the at¬ 
tenuation curve fou nd tor local UV starburst galaxies 
( CaUetti et al |2000 ), 1 )ut its presence to varyi ng degrees 
has been noted in higher-redshift galaxies (e.g., Noll et al. 


20091 Buat et al.|2011 2012 Kriek fc ConroyJ26l3HScov~ 

ille et al.||2015|). 

'id ascertain the existence of this excess absorption, we 
examined the average photometry of galaxies in bins of Tb 
and E(B—V), where the latter were dete rmined from the 
SED modeling discussed in Section [T6| To improve the 
rest-frame wavelength sampling of the average photom¬ 
etry, we include all galaxies irrespective of their sSFRs. 
We then calculated the average, and error in the average, 


flux densities in rest-frame wavelength bins of 100 A(e.g., 
similar to the method used above to compare the average 
photometry of galaxies in our sample to the functional 
fit to the photometry; Figure]^. A linear function was 
fit to the log of the average flux density as a function of 
log(A) in two wavelength ranges bracketing the 2175 A 
feature: A = 1500 — 1800 and 2600 — 3000 A. The errors 
in the mean flux densities were used to weight the points 
in the fitting. The average flux densities throughout the 
wavelength range from 1500 to 3000 A were then divided 
by this linear fu ncti on to produce the normalized points 
shown in Figure 

There is on average a lower flux at wavelengths corre¬ 
sponding to roughly where the excess absorption would 
be expectedp^ While the depth of the absorption ap- 


As the wavelength windows used to compute the UV slope 
largely exclude the region where we observe a deficit in the average 


pears to be insensitive to Tb, there is some marginal sig¬ 
nificance for a trend between the absorption depth and 
the color excess of the stellar continuum. A Spearman 
correlation test implies that that the null hypothesis of 
no correlation between the depth of the 2175 A and the 
color excess of the stellar continuum can be ruled out at 
the ~ 3cr level. We caution that the excess absorption 
will be smoothed by the finite wavelength sampling of the 
average photometry as well as the fact that the average 
photometry is computed from broad- and intermediate- 

band magnitudes. _ 

For comparison, Buat et al. (2012) deduce the shape of 
the dust attenuation curve for a sample of 751 galaxies at 
0.95 < 2 < 2.2, and find evidence for excess absorption at 
2175 A for roughly 20% of their sample, almost all (90%) 
of which lie at z < 1.5. They also find that the rate of 
a secure detection of bump increases by a factor of « 2 
when considering the dustier galaxies (i.e., those directly 
detected at mid- and far-IR wavelen gths). In another 
recent study, |Kriek & Conroy (2013) examine a sample 
of somewhat lower-redshitt (0.5 < z < 2.0) galaxies and 
find that the strength of the 2175 A absorption correlates 
with the “steepness” (i.e., the wavelength dependence) of 
the attenuation curve. However, we note that the galax¬ 
ies in their sample have Wua < 140 A, values that are 
substantially lower (corresponding to lower sSFRs) than 


those considered here (Whq > 100 A). More definitive re¬ 
sults on the extent to which the 2175 A absorption may 
be present in the galaxies targeted with MOSDEF, as 
well as the characteristics (e.g., slope of the attenuation 
curve, sSFR, redshift) upon which its depth may depend, 
will come from the larger (completed) survey dataset. 

Regardless, our analysis shows that the attenuation 


normalized flux, the effect of imposing a linear fit to the flux as 
a function of wavelength will result in a minimal bias of A0 < 
0.1. This small bias does not affect our conclusions regard ing th e 
contribution of A stars to the near-UV continuum (Section|3.3.2|l, 
though of course fitting the slope without the photometry that may 
be affected by the 2175 A absorption will result in larger random 
uncertainties. 
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Figure 12. Normalized average photometry as a function of rest- 
frame wavelength for galaxies in bins of Balmer optical depth 
(top) and color excess of the stellar continuum (bottom). The av¬ 
erage photometry of galaxies is computed in bins of rest-frame 
wavelength of AA = 100 A, and the average photometry is nor¬ 
malized using a linear fit to the points over wavelength ranges 
A = 1500 - 1800 and 2600 - 3000 A. 

curve for the MOSDE F galaxies exhibits a slope that is 


identical to that of the Calzetti et al. ( 2000 ) attenuation 
curve for A < 0.6/im, but with a lower normalization 
(i?y). For A > 0.6/im, the MOSDEF curve is similar 
to that of the SMC curve, and more generally the two 
are essentially identical for A > 0.25 ^m. Recent stud¬ 
ies have suggested that some populations of high-redshift 


Reddy et al. 1120061 Siana et al. 

120081120091 

Reddy et al.j 

2010||Shim et al.|2011 Oesch et al.|2013). h 

lere, we have 


shown that through direct measurements of the attenua¬ 
tion curve, an SMC-like behavior at longer wavelengths 
(A > 2500 A) may be more generally applicable at high 
redshift. A comparison between physical models of the 
distribution of dust and stars in galaxies and the attenua¬ 
tion curve will be presented elsewhere, but we note that 
the shape of the attenuation curve may imply a larger 
covering fraction of dust in the rest-optical, or a differ¬ 
ent du st grain size distribut ion, than would be indicated 
by the Calzetti et al.|(2000l curve. The implications for 
thi 


these ditterences in the shapes of the attenuation curves 
are discussed in the next section. 

4. DISCUSSION AND IMPLICATIONS 
4.1. Reassessment of the Stellar Populations 

Given the apparent differences in the attenuation curve 
derived in this work relative to the cur ve typically as¬ 


sume d for high-redshift galaxies (i.e., Calzetti et al. 
20001 , we adopted our new calculation of the attenuation 


curve to refit the stellar populations of galaxi es in our 
sample, in the same way as described in Section [2?6l The 
comparison of several stellar population param eters us¬ 
ing our new curve versus those obtained with the Calzetti] 
et al. (2000) and SMC curves is shown in Figure 13 The 
mean color excess of the stellar continuum assuming the 
MOSDEF attenuation curve is (A E(B — V)) « 0.02 red¬ 


der than that obtained with the Calzetti et al.j (2000| 
relation {E{B — V) « 0.19 vs. 0.17). Combining this av¬ 
erage difference in the color excesses with the Afc « 1.63 
difference in the two attenuation curves at A = 1600 A we 
infer dust obscuration factors (required to de-redden the 
observed UV SFRs) that are « 1.2 times smaller 


than those inferred from the Calzetti et al. (2000) rela¬ 


tion. As a result, the average SED-inferred SFRs ob¬ 
tained with the MOSDEF attenua tion curve are » 20% 
smaller than those obtained with a Calzetti et al. (2000| 
curve. This systematic difference in the Sh'Rs is typically 
smaller than the random uncertainties, where the latter 
were computed using Monte Carlo realizations of the ob¬ 
served photometry and refittin g the model SEDs to the 
perturbed photometry; Section 2.6). Similarly, the sys¬ 
tematic offset seen here is well within the typical factor 
of « 2 uncertainty in independent measurements of the 
SFRs based on combining UV and IR data for galaxies 


redshifts, as those considered here ( 

Reddy et al. 

2006 

Daddi et al. 

2007] Reddy et al. 

2010 

2012a|). These re- 


the difference in E{B—V) implied by the two attenuation 
curves is a function of the “redness” of the object. For 
the bluest galaxies, we find AE{B — U) « 0 (these are 
galaxies for which little reddening is required to repro¬ 
duce the UV colors). For the reddest galaxies (e.g., those 
with E{B — V) > 0.3 as computed from the MOSDEF 
attenuation curve), we find AE{B — U) « 0.08. 

Assuming the MOSDEF attenuation curve results in 
galaxies that are on averag e » 300 Myr older rela tive to 
the ages derived using the jCalzetti et al.j (|2000|) curve, 
but we view these differences to be negligible given the 
large uncertainties in the agesj^ Note that the “uncer¬ 
tainties” quoted here do not account for other system¬ 
atic effects, including differences in the assumed star- 
formation histories of galaxies. Finally, we note that as¬ 
suming the MOSDEF attenuation curve yields better fits 
to the longer wavele ngth (i.e., rest-frame ~ 1.0 — 2.5 /im) 
photometry than the Calzetti et al (2000| curve for rising 
star-format ion histories (1'igure T?). As such, the stellar 
masses are « 25% lower whe n we adopt the MOSD EF at¬ 
tenuation curve in lieu of the Calzetti et al. (2000) curve. 
This implies an « 0.1 dex difference in log ill that is of 
the same order as the uncertainty in stellar mass implied 
by the photometric errors. Similarly, in comparison to 
the results derived with an SMC attenuation curve, those 
assuming the MOSDEF curve imply significantly redder 
E(B-V) {AE{B — V)~ 0.10), younger ages (by a factor 
of « 2), higher SFRs (by a factor of « 2), and lower stel¬ 
lar masses {A\og{M*/M q) « 0.16dex). The differences 
stem from the fact that the SMC curve has a steeper 
wavelen gth dependence i n the UV than either the MOS¬ 


DEF or Calzetti et al. (2000) attenuation curves, and 
thus a lower amount of reddening is required to fit a given 
UV color. As a result, the SFRs are lower with the SMC 
curve. Further, the near-IR-optical color is attributed 
more to an older stellar population than dust reddening, 
resulting in older ages and larger stellar masses. 

Another point of consideration is how the sSFRs are af¬ 
fected by different assumptions of the attenuation curve. 
Recall that the sSFRs were calculated by combining Ha- 
based SFRs, which were determined independently of 


A subset of objects have “best-fit” ages that are a factor of 
> 2 older with th e MOSDEF attenuati on curve relative to those 
obtained with thelCalzetti et al.l ll2000t curve. T hese are objects 
where the MOSDUh and |(Jalzetti et al.| (|2000|l “best-fit” SEDs 
are essentially identical in tne rest-lrame optical, but where the 
former are generally fainter in the rest-frame near-IR. As such, 
the implied Balmer/4000A breaks are stronger for the MOSDEF 
“best-fit” SEDs resulting in older ages. 
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Figure 13. Comparison of the color e xcesses, ages, SFRs. a nd stellar masses returned from SED fitting assuming our new attenuation 
curve relative to those obtained with the|Calzetti et al.|(|2000|l and SMC attenuation curves. The solid lines denote a one-to-one agreement 
between the plotted quantities. 
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Figure 14. Residuals (i.e., difference in the observed and model 
photometry) for galaxies in our sample assumin g a rising star- 
format ion history and the MOSDEF {cyan) and jCalzetti et al.| 
l|2000|| (orange) attenuation curves. The mean residuals and the 
errors in these means are in dicated by the blue a nd red symbols 
assuming the MOSDEF and jCalzetti et al.| | |200^ curves, respec¬ 
tively. For rest-frame wavelengths between ps i and 2.5/i.m, The 
MOSDEF attenuation curve results in a statistically significant im¬ 
proved fit to the obser ved photometry , on av erage, relative to that 
obtained assuming the jCalzetti et ^ | |200l5| l curve. 


SE P fit ting, with the SED-based stellar masses (Sec- 
tion|2,6 ), The Ha-ba sed SFRs were derived assimhng the 
1989) extinction curve (Section [2^^ , The 


Cardelli et al. 


SEP-based ste 


lar ih asses were computed assuming the 


Calzetti et al, (2000) curve. We will return to a more de- 
tailed discussion ot the i mplic ations of our results for SFR 
determinations (Section 4.4), but for the remainder of 
our analysis, we adopt the sSERs computed assuming the 
stellar masses obtained with the MOSPEF attenuation 


curve, and SFR(Ha) obtained with the Cardelli et al. 
(1989) extinction curve. As shown in Figure TS the new 
sSFRs computed in this way a re quite similar to th ose 
computed when we assume the Calzetti et al. (2000) at¬ 
tenuation curve in deriving the stellar masses, with the 
former being Alog(sSFR/yr“^) ~ 0.09dex smaller than 
the latter. Because there is only a systematic offset be¬ 


tween the two sSFR determinations, and the relative or¬ 
dering of the galaxies is preserved, the sS FR b ins used 
to calculate the attenuation curve (Section 3.6.1) will in¬ 
clude essentially the same set of galaxies regardless of 
which attenuation cu rve is assumed (and, in any case, 
the results of Section |3.6.1| suggest that because the at- 
uted 


tenuation curves computed from the two sSFR bins are 
very similar, we could have simply ignored any binning 
with sSFR in computing the curve). 



sSFR (yr q - M* from Calzetti 


Figure 15. Specific SFRs calculated assuming Ha-based SFRs, 
and stellar masses from SED fitting. SFR(Ha) assumes the 
iCar delli et al. (|1989l> ext inction curve. The stellar masses assume 
the|Ualzetti et al.| (|2000|l (ordinate) and the MOSDEF (abscissa) 
attenuation curves. 


4.2. Measurements of UV Slopes 

/3phot(1260 — 1750 A) is not 


As noted in Section 


3.3.2 


substantially bluer than /3phot(1260 —2600 A), suggesting 
that A stars do not contribute significantly to the UV 
continuum for the galaxies in our sample. It is worth¬ 
while to consider this point further, because an intrinsi¬ 
cally blue spectral slope in the range 1260 — 1750 A can 
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appear redder than the slope measured over the wider 
baseline in wavelength (1260 —2600 A) if the attenuation 
curve is rising steeply in the UV. Therefore, we must rule 
out the case where the attenuation curve masks the con¬ 
tribution from A stars to the near-UV continuum. Our 
new derivation of the attenuation curve implies that the 
difference between /3(1260—1750 A) and /1(1260—2600 A) 
for galaxies where the latter is < 0.5 (a limit that en¬ 
compasses almost all of the galaxies in our sample), is 
A/3 < 0.2 for E{B—V) < 0.4. This difference in UV slope 
is small compared to both the typical measurement errors 
in the UV slopes and the dispersion in /3(1260 — 1750 A) 
at a given /3(1260 — 2600A) (Figure]^. Thus, the red¬ 
dening implied by the MOSDEF attenuation curve is in¬ 
sufficient to make an intrinsically blue ,3(1260 — 1750 A) 
appear substantially redder than /3(1260 — 2600 A), and 
our conclusion regarding the contribution of A stars to 
the UV continuum still holds. 

The test just discussed allows us to calculate the useful 
relationship between /3 and the color excess of the stel¬ 
lar continuum assuming the MOSDEF attenuation curve 
and an underlying solar metallicity stellar population 
with constant star formation for at least 100 Myr (i.e., 
with an intrinsic spectral slope of /3int ~ —2.44): B = 
-2.44-1-4.54 X E(B — V). Additionally, Equationsand 
[pimply that the attenuation in magnitudes at 1600 A is 
related to the color excess by = 8.34 x E(B — V). 

Thus, combining the previous two equations, we find 


A * = 1.84^ -f 4.48, 
1600 A ^ ’ 


(9) 


again for an intrinsic spectral slope /3int — — 2 .44. F or 


comparison, the relation from Meurer et al. (19991 is 


nen 


A s = 1.99/3 -I- 4.43, or A s = 1.99/3 -I- 4.86 w 
1600 A ’ 1600 A ^ 

shifted to account for the difference in the intrinsic UV 


slope (/3o) between Meurer et al. (1999) and the present 
study (i.e., /Bg = —2.T6 vs —2.44). Consequently, for the 
bulk of the galaxies in our sample with /3 < —0.5, the im¬ 
plied difference in the magnitudes o f extinction at 1600 A , 
when assuming the MOSDEF and Meurer et al. (19991 
attenuation curves, is AA s < 0.31. More specih- 
cally, the MOSDEF attenuation curve implies dust ob¬ 
scuration factors, 10"[0.4xA(1600 A)], that are a factor 
of 5. 1-32 tim es lower than those imp lied by the Meurer 
et al.| (|1999|) (or jCalzetti et al.||2000 ) attenuation curve, 
with an average correction that results in SFRs that are 
« 20% lower for the MOSDEF curve, as noted above. 
Equation!^ can be used to directly determine the attenu¬ 
ation of the UV continuum for galaxies of a given ( mea- 
sured) UV slope. However, as we discuss in Section |4.4[ 
additional factors may be required to recover total SFHs 
from UV-based ones. 


4.3. Comparison of the Color Excesses and Total 
Attenuation of the Stellar Continuum and Ionized 
Gas 

There has been a considerable effort to investigate the 
relationship beteween the attenuation of the stellar con¬ 
tinuum and ionized gas in high-redshift galaxies, par¬ 
ticularly in li ght of local observations of star-forming 


galaxies (e.g., Fanelli et al.| 19881 Calzetti et al. 

1994 

|Mas-Hesse &: lvunth|1999 Moustakas & Kennicutt 



Kreckel et al. 2013) that suggest a higher reddening of 
the ionized gas than of the stellar con tinuum : e.g. 


EjB — U)stars = 0.44 X E{B — U)gas (Calzetti et al. 


2000). Note that this often-quoted relationship describes 


the color excess of the stellar continuum and the nebular 
line-emitting regions, not the total attenuation. More 
precisely, this relation is derived under the assumption 
that the starburst attenuation curve applies to the stel¬ 


lar continu um, and a Galactic extinction curv e (Cardelli 
et al.||198'9 ) to the ionized gas ( Calzetti||2001 ). However 


most investigations of the stellar and nebular reddening 
in high-redshift galaxies have assumed that the same at¬ 
tenuation curve applies for both the gas and the stars, 
and it is unclear if this is an appropriate assumption in 
the context of high-redshift galaxies. 

From the standpoint of local galaxies, the difference in 
attenuation curves can be understood in the framework 
where the recombination emission arises from spatially 
compact HII regions that are subject to a high cover¬ 
ing fractio n of dust and where a line-of-sight extinction 


curve (e.g., Cardelli et al. 1989 or SMC) is more appropri¬ 
ate, versus the more spatially extended stellar continuum 
that has a lower dust covering frac tion and is subject t o 


a “greyer” attenuation curve (e.g., Calzetti et al. 2000). 


This conf iguration is simi l ar to t he two component dust 


model of Chariot & Fall (2000) where the HII regions 
are attenuated by a component of dust in undissipated 
parent birth clouds, with a second diffuse component of 
dust in the ISM. 

Studies of UV-selected galaxies at z ~ 2 have shown 
a general agreement between UV and Ha-based SFRs 
when applying the same color excess and assuming the 


Sas (e.g., 

Erb et al.| 2006a 

Reddy et al. 

2010 

Cowie & Barger 


tor results at z < 0.5). 


highly star-forming or more massive galaxies, and/or 
photometrically-selected galaxies, indicate a higher at¬ 
tenuation of the gas versus the stellar continuum, again 
with most studies assuming that the same atte nuation 


Schreiber et al.|2009| Yoshikawa et al. 

20101 IWuyts et al. 

2011 

IKashino et al.|2013[ |Wuyts et a 

'.||2013t|Price et al. 


2014). 


4.3.1. Comparison of Color Excesses 

A plausible scenario that can resolve the current obser¬ 
vations is one in which the difference in the continuum 
and nebula r attenuation is a function of SFR, m ass, or 
sSFR (e.g., Reddy et al.J|2010 Price et S(]|2014 ). As re¬ 
marked in Section |3. 6. 2f the selective attenuation curve 
derived in this study implies that on average the ionized 
gas suffers greater extinction than the continuum (at the 
wavelengths of Ha and H/3). We scrutinized this result 
by comparing the color excesses of the stellar continuum 
derived from SED fitting, with those of the ionized gas, 
computed as: 


E{B - U)gas = 


2.5 

fc(H/3) — fc(Ha) 
1.086 

fc(H,0) — fc(Ha) 


/Ha/H/3\ 
2.86 ) 

Tb ~ 0.950rb. 


( 10 ) 
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Figure 16. Left: Compari s on of the color excesses derived for the stellar continuum and the ionized gas. E(B — F)gas is computed 
assuming the|Cardelli et al.|||1989|l extinction curve, and E{B — F)stars is the value returned from the SED fitting when we assume the 
MOSDEF attenuation curve. Ihe solid line denotes identical colo r excesses derived for the stellar and gas components, and the dashed 
line indicates the relation E(B — F)gas = E(B — F)stars/0.44 from |Calzetti et al.| (|2000|l. Points are differentiated according to the sSFR. 
The large star denotes the average values for the H/3-undetected galaxies. Higtit: Uitterence between gas and continuum color excesses as 
a function of sSFR, where the solid line indicates no difference between the color excess of the nebular regions and the stellar continuum. 
The large star denotes the average values for the H/3-undetected galaxies. 


In the last step, we have assumed that the Cardelli 


et al. (1989) extincti on curve applies to the HII re- 


gions. Adopting the Calzetti et al. (20001 and MOS¬ 
DEF attenuation curves results in — V)gas values 
tha t are a factor o f 0.90 x and 1.01 x those assuming 
the iCardelli et al.| (|1989|) curve, respectively. There is 
a large scatter m the relationship between the gas and 
continuum color excesses, with « 50% of the objects hav¬ 
ing equal gas and continuum color excesses within the 
(1 a) uncertainties, and only 7% consistent with having 
E{B — E)gas > E{B — E)stars/0.44. As could have been 
anticipated from Figures and the degree to which 
E{B — F)gas diverges from E{B — F)stars is dependent 
upon the sSFR of the galaxy, with the largest differences 
for galaxies with the largest sSFRs (Figure 16) ^ 

We should consider, however, whether sel^ion biases 
could result in correlation between the difference in the 
gas and nebular color excesses and the sSFR. In particu¬ 
lar, if galaxies with Ha fluxes close to our detection limit 
of ~ 5 X 10“^® ergs“^ cm”^ are dusty, then their H/3 lines 
may be undetected with a greater frequency than galax¬ 
ies of the same dust content but with bright Ha lines. As 
a result, we may be insensitive to galaxies with low SFRs 
and high nebular color excesses. However, as seen in Fig¬ 
ured the H/3-undetected galaxies in our sample have lim¬ 
its that span the full range of observed Lhq- Further, if 
we compute the average SFR(Ha), sSFR, E{B — H)gas, 
and E{B — Hjst. ar.s fo r the H/3-undetected galaxies, as 
shown in Figures 16 andl^ (where the averages for the 
SFR and E{B — V )gas areaetermined from the compos¬ 
ite stack of their spectra; Section]^, the aforementioned 
quantities suggest that there is not a large population 
of dusty galaxies with low SFRs and low sSFRs. Aside 
from this direct observation, it has been shown in sev¬ 
eral studies that the mean obscur ation of high-redshift 


galax i es increases with SFR (e.g., Adelberger fc Steidel 


2000[ [Reddy et al. 2006 Huat et "aT. ' 2007 ~ Hurgarella 


As we sh ow in Sectfion |4.4.2[ a different ^lati on is found 

between E(B — P jgas — E(B — V )sta.rs and s8FR(8FU), i.e., the 

sbf K compufed wtien we use SEU^Uased Bl* 


et al. 

2009 

Reddy et al. 

2008 

Buat et al. 2009 Reddy 

et al. 

p610 

Sobral et al. 

2012 

Dominguez et al.||2013 

see also discussion below), suggesting that among galax- 


ies with low Ha fluxes, those that are intrinsically less 
dusty vastly outnumber those with high dust content. 
The stacking of mid- to far-IR data for the MOSDEF 
galaxies will help to clarify these arguments. There is 
no similar bias in terms of E{B — H)stars) as the spec¬ 
troscopically confirmed MOSDEF galaxies have rest-UV 
colors similar to t hose of galaxies in the parent sample 
( Kriek et al.||2014 ). 


Formally, a Kendall’s r test indicates a probability of 
< 1.8 X 10“^ that the difference between gas and stellar 
color excesses is uncorrelated with sSFR. The correla¬ 
tion with Ha-based SFR is even stronger, with < 10“® 
probability that E{B — VXas — E{B — H)stars is uncorre¬ 


lated with SFR (Figuredl7|). Similarly, there is a < 10“"^ 
probability that the difference in color excesses is uncor¬ 
related with stellar mass. Thus, the conflicting results 
on the relationship between the color excesses of the HII 
regions and stellar continuum in high-redshift galaxies 
can be reconciled if selection over different ranges of SFR 
and/or stellar mass leads to the inclusion of more or fewer 
galaxies with the largest discrepancies between nebular 
and continuum color excess. Specifically, galaxies with 
SFRs> 20MQyr-^ exhibit nebular color excesses that 
are substantially redder than the continuum color ex¬ 
cesses. 

In contrast, over the range of SFRs characteristic of 
UV-selected galaxies at z ~ 2 (< 2Q for a 


Chabrier IMF; e.g., Reddy et al. 2012b), the nebu¬ 


lar color excesses are essentially equivalent to or bluer 
than those of the stellar continuum. Of the 51 ob¬ 
jects that have E{B — K)gas that is more than 1 cr 
bluer than E{B — K)stars, 35 (or « 69%) formally have 
E{B — V)gas = 0. These galaxies are unlikely to have 
HII regions that are compl etely dust-free given the rea¬ 
soning presented in Section 3.5 and their E(B — H)stars 
span the full range observed for galaxies in our sample 
in general. An analysis of the spatially-resolved colors of 
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these galaxies should lend more insight into the question 
of why they exhibit E{B — y)gas that are bluer than the 
color excess of the stellar continuum. In summary, we 
find that the difference in the color excesses of the stellar 
continuum and ionized gas in high-redshift galaxies is a 
strong function of SFR: above an SFR of « 20 Mq yr“^, 
the nebular lines become increasingly reddened relative 
to the color excess of the stellar continuum. 

4.3.2. Comparison of the Total Attenuation of the Ionized 
Gas and Stellar Continuum 

So far, our discussion has focused on a comparison of 
the color excesses of the ionized gas and the stellar con¬ 
tinuum in high-redshift galaxies. Recall that the color 
excess of the i onized gas was computed assuming the 
Cardelli et al. (1989) extinction curve. The color excess 
ot the stellar continuum was derived assuming the MOS- 
DFF attenuation curve. These respective curves were 
used to also compute the total attenuation in magnitudes 
of the line and continuum emission at 6565 A: 

Ahq = fccardem(6565A) x£;(R-F) gas 

= 2.52xF(R-F)gas; 

Acont = fcMOSDEF(6565 A) X E{B — V) stars 

= 1.92xF;(R-F)stars. (11) 


The comparison of Ana an d Acon t (Figure 181 confirms 
the result noted in Section |3.6.2[ namely mat the ion¬ 
ized gas is on average more obscured than the stellar 
continuum. Specifically, we find that the average dif¬ 
ference in the magnitudes of attenuation between re¬ 
combination line and continuum photons at 6565 A is 
(Anct — Acont) — 0.15. However, as was the case in the 
comparison of the color excesses, we find that the differ¬ 
ence in total attenuation of line and continuum photons 
depends on the SFR, such that galaxies with the largest 
SFRs can exhibit line-emitting regions that are half a 
magnitude (or more) attenuated than the stellar contin¬ 
uum. 

4.3.3. Physical Interpretation 

A final point of consideration is the stron g de pendence 
of total dust attenuation on SFR (Figure 20). Such a 
correlation has been noted before both awow (Wang 
fc Heckm an||1996| [Hopkins et al.||2001|) and high red- 
shift (e.g.,|Adelberger & Steidel|2U0Ut IRe 


2008 


2012 


ger fc Steidel|20U0| [Reddy et al.|2006 

Buat et al.||20U7[ [Hurgarella et al.||20U9[ [Reddy et al. 


IBuat et aT 12009 [Reddy et al.||2010 [Sobral et al. 

Dominguez et al. 2013[), and has alsobeen shown to 


evolve~with redshitt (e.g., Reddy et al.||2010 Dominguez 


et al. 2013). Physically, this trend can be understood 
as the product of increased metal and dust enrichment 
of the ISM with increasing SFR. In the context of the 
present study, a possible simple model that can account 
for the trends in ionized gas versus continuum attenua¬ 
tion is one in which a modestly reddened stellar popula¬ 
tion dominates the observed UV through optical contin¬ 
uum emission of galaxies over the entire range of SFRs 
considered here. Meanwhile, as the total SFR rises, a 
greater fraction of this SFR is obscured in optically-thick 
regions, as a consequence of the increased dust enrich¬ 
ment that accompanies galaxies with larger SFRs. As a 
result, the average optical depth towards ionizing stars 


rises with SFR. This simple picture is illustrated in Fig- 
We emphasize that this is an “average” effect, as 
ear that some fraction of O stars are in relatively 
unobscured regions of the galaxies where they can con¬ 
tribute significantly to the UV continuum and where the 
average reddening (as measured by the UV slope, j3) is 
lower than that measured from the Balmer lines. 

The framework just discussed is somewhat different 
than the one originally invoked to explain the dispar¬ 
ity between the gas and continuum color excesses in 
the local universe. Locally, such a discrepancy between 
color excesses is thought to arise from the increased dust 
covering fraction surrounding HII regions (e.g., perhaps 
due to undissip ated parent birth clo uds surrounding the 
youngest stars; Calzetti et al. 1994). If this is the case 
at high redshift, i.e., that most ot the short-lived O stars 
reside in or close to their parent molecular clouds, then 
we should see a systematic offset between gas and contin¬ 
uum color excesses at all SFRs. However, this expected 
systematic offset is not seen in our data (e.g.. Figures 16 


17 and 18). Rather, our results suggest two distinct steL 
lar populations in each galaxy: one in which stars of all 
masses are subject to modest amounts of dust redden¬ 
ing and which dominate the observed UV through op¬ 
tical continuum luminosity; and another in which stars 
of all masses reside in very dusty regions. The second 
population becomes more dominant over the first pop¬ 
ulation (in terms of its contribution to the nebular line 
and bolometric luminosities) as the total SFR of a galaxy 
increases. Consequently, the color excesses (and total at¬ 
tenuation) of the gas and continuum diverge with rising 
SFR. This simple picture may also imply that the atten¬ 
uation curve relevant for the ionized gas may transition 
from a gray curve (e.g., Calzetti et al. 2000) to a pro¬ 
gressively “steeper” (e.g.. Galactic or SJVIC) curve with 
increasing SFR. 

Along these lines, more recent high spatial resolu¬ 
tion observations of several nearby galaxies suggest that 
the relationship between the attenuation of the gas and 
stars varies with SFR surface density (X)sfr)i such that 
the Balmer lines are systematically more attenuated 
than the stellar continuum for regions w here ^ 

O.OI M c^ yr“^ kpc~^ (e.g., see Figu re 4 of Kreckel et al. 
2013 see also Boquien et al.||2015). These observations 
have been interpreted as evidence that the Balmer lines 
are dominated by du st-buried HII regions at high X)sfr 
(K reckel et al. 2013). If the higher SFR galaxies in our 
sample also have higher X)sfr’ then the situation in¬ 
ferred locally on small scales (e.g., 100-200pc) may also 
apply at high redshift. This issue can be addressed with 
the larger MOSDFF sample, which will include galaxy 
size information. 

Our analysis demonstrates that while the wavelength 
dependence of the dust obscuration appears to be in¬ 
dependent of specific SFR, the total obscuration of the 
ionized gas relative to that of the stellar continuum 
varies with SFR and sSFR. Specifically, galaxies with 
larger SFRs (and sSFRs) exhibit redder color excesses 
and larger total attenuations for the line emission rela¬ 
tive to those for the continuum. These trends are re¬ 
sponsible for the scatter seen between /3 and (Fig¬ 
ure [^. The physical picture discussed here implies that 
galaxies with larger SFRs may exhibit a higher disper- 
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Figure 17. Difference between gas and continuum color excesses as a function of Ho SFR (left) and stellar mass (right), where the solid 
lines indicate no difference between the color excess of the nebula r regions and the stell ar continuum, and the large stars indicate the values 
for H/3-undetected galaxies. SFR(Ho) is computed assuming the |Cardelli et ahj l llQS^ extinction curve. 



Figure 18. Difference between the magnitudes of attenuation of 
recombination line photons and continuum p hotons at 6565 A. as a 
function of SFR. Ahq a,nd Acont assume the|Cardelli et al.|(|1989ll 
Galactic extinction curve and the MOSDEF attenuation curve, re¬ 
spectively. The solid and dashed lines denote equal attenuation 
and the mean difference ((Aho ~ Acont) — 0.15) in attenuation of 
line and continuum photons, respectively. 

sion in their spatially resolved colors as the dust obscu¬ 
ration becomes more patchy with the rise of a second 
(dustier) population. Alternatively, the increased dust 
covering fraction expected for galaxies with larger SFRs 
may be reflective of a compact configuration of the star 
formation. Hence, future prospects for elaborating on 
this simple picture will include a detailed analysis of 
how the SFR surface density, and the spatially resolved 
UV colors of galaxies—whose measurement is aided with 
deep high-resolution HST imaging in the fields targeted 
with MOSDEF—correlate with globally-measured color 
excesses of the ionized gas and the stellar continuum. 


4.4. Implications for SFR Measurements 
4.4.1. Ha and UV-based SFRs 


As mentioned in Section |4.1[ the SED-based SFRs 
assuming a MOSDEF attenuation curve are on aver- 
age lower (by » 20% ) than those determined with the 
Calzetti et al.| (|2000| curve. The Ha-based SFRs also 
rely on an assumption of the attenuation curve, both in 
deriving the color excess of the ionized gas, and in apply¬ 
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Figure 19. Illustration of a simple geometry of dust and gas 
that can account for the trends in the difference between ion¬ 
ized gas and continuum color excess, and total attuation, ver¬ 
sus SFR. The yellow region denotes the diffuse dust component 
(that may be patchy). The red regions indicate areas of increased 
dustiness within the galaxy. The blue and red stars indicate high 
mass (ionizing) and lower mass stars, respectively. At lower SFRs 
(< 20 Mq yr~^), stars of all masses are uniformly obscured. As the 
SFR increases, the diffuse component becomes more dust-enriched 
(as indicated by the darker shade of yellow), while regions of more 
highly obscured SFR (red regions) become prominent. As the SFR 
increases, these more obscured regions begin to dominate the neb¬ 
ular line and bolometric luminosities. The diffuse component dom¬ 
inates the UV through optical SED at both low and high SFRs. 


ing the total dust correction to recover the intrinsic Ha 
SFRs. The d efault values o f SFR (Ha) were computed 
assuming the Cardelli et al. ([1989) extinction curve — 
i.e., this curve was used to calculate E{B — V)gas given 
the Balmer decrement, and the value of the curve at 
6565 A was used to dust correct the observed Ha lu- 


minosity (Section 

2.8 

). The more common practice of 

assuming that the 

Ca 

Izetti et al. 

(2000 

) curve applies to 

high-redshilt Hii regions results 

n Sb'Ks that can be as 


much as » 5 0% larger than those assuming the Cardelli 
et al. (1989) extinction curve (Figure [21]), with an aver¬ 
age offset of « 20%. This comparison underscores the 
importance of considering the appropriate attenuation 
curve to use when dust correcting the line and contin¬ 
uum emission in high-redshift galaxies. 

When compared against each other, the Ha and SED- 
determined SFRs show a ge neral agreement, with a scat¬ 
ter of 0.34 dex (Figure [2^. A sim ilar agreement has 
been noted in Steidel et al. (2014), who also compare 
Ha and SED-based SFRs, but who use a different pre¬ 
scription for estimating the color excess relevant for the 
nebular emission as they did not have Balmer decrement 
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Ha/H/3 

Figure 20. Ho-based SFR versus Balmer decrement. The large 
star denotes the average values for the H/3-undetected galaxies. 
Note that the Balmer decrement is used in the calculation of the 
SFR, resulting in a tightening of the scatter shown here. Nonethe¬ 
less, the results suggest the same trend between SFR and dustiness 
that has been noted in other studies at low and high redshift (see 
text). 



Figur e 21. Ratio of Ha-based S FRs assuming t he ICalzetti et al.| 
(|2000|l attenuation curve and thelCa rdelli et aril|l989ll ex tinction 
curve, versus the those assuming IhelUardelli et aT] (|iy89|| extinc¬ 
tion curve. The large star denotes the a verage values tor the H/3- 
undetected galaxies. Assuming that the|Calzetti et al.|(|2000|l at¬ 
tenuation curve applies to the HII regions ot tiigti-reclshitt galaxies 
results in dust-corrected SFR(Hq;) that can be as much as ^ 50% 
larger than those derived assuming the |Cardelli et al.| ||1989|l ex¬ 
tinction curve. 

measurements for their sample. Notwithstanding this 
general agreement, the fraction of objects for which the 
Ha-based SFR exceeds the UV- (or SED-) based SFR 
by more than 1 a increases from / = 0 to / = 0.8 
from SFR(Ha)~ lM 0 yr“^ to 3 OOM 0 yr“^. We for¬ 
malized this trend by fitting a third-order polynomial to 
SFR(Ha)/SFR(UV) versus x = log[SFR(Ha)/A /0 yr“^]: 

= 0.29 + 0.31a; + 0.293;^ + 0.0040x^ (12) 

SFR(UV) 

Figure]^ shows that > 40% of galaxies with SFR(Ha)> 
2 OM 0 yr“^ have UV- and SFD-based SFRs that under¬ 
predict the Ha-based SFR by an average factor of ~ 2. 
It is unlikely that this trend is due to selection effects, 
as the MOSDFF galaxy selection is based primarily on 


rest-optical continuum emission, and the success of iden¬ 
tifying a redshift depends only on whether any of the 
rest-optical emission lines are detected, independent of 
the UV continuum. Thus, it is improbable that we are 
missing a large number of galaxies where the UV-based 
SFR is substantially smaller than the Ha-based SFRs for 
SFR(Ha)< 20Moyr-h 

Focusing on those galaxies with lower SFRs, we find 
that their SFR(Ha)/SFR(UV) ratios lie systematically 
below a value of unity. We must address whether this 
systematic difference is due to variations in the dust 
attenuation of the nebular regions and stellar contin¬ 
uum, or if there is some other mechanism that may be 
partly responsible for the low SFR(Ha)/SFR(UV) ra¬ 
tios for these galaxies. To investigate this, we calculated 
the ratio of the observed Ha and UV SFRs (i.e., un¬ 
corrected for dust extinction) for all the objects in our 
sample (Figure 1^). The observed SFRs will include 
the light emitteoby unobscured stars, as well as the 
unattenuated light from obscured stars. Galaxies with 
(dust corrected) SFR(Ha)< 2 OM 0 yr“^ have observed 
SFR(Ha)/SFR(UV) ratios that are generally within 3 a 
of a ratio of unity. As a consequence, for the stars con¬ 
tributing to the observed Ha and UV luminosities, there 
is no evidence that we are observing these them several 
tens of Myr after an initial burst of star formation: this 
scenario would result in observed SFR ratios that devi¬ 
ate from unity due to the different timescales (and mass 
ranges) over which the Ha and UV luminosities are sen¬ 
sitive to the total SFR. Thus, if the stars dominating 
the unobscured Ha and UV luminosities have a similar 
IMF and star-formation history as those stars that dom¬ 
inate the obscur ed Ha and UV luminosities, then the 
results of Figure suggest that the low dust-corrected 
SFR(Ha)/SFR(UV) ratios measured for galaxies with 
SFR(Ha)< lOM 0 yr“^ are likely tied to differences in 
the dust obscuration of the nebular emission and stellar 
continuum, rather than differences in star-formation his- 
tory, changes in the IMF, or sampling stochasticity (c.f., 
Lee et al.|2009[|Meurer et al.|2009[[Fumagalli et al.|20Il[ 


!talJ2009l|Meurer i 

ilva et al. [2014[ Ha 


ommguez et al. 112014 ( 


4.4.2. A Prescription for Estimating E{B — U)gas and Ahc 


With our independent measurements of the color ex¬ 
cesses of the stellar continuum and ionized gas, we can 
quantify the relationship between SFD-inferred quanti¬ 
ties (i.e., E{B — U)stars and SFR(SFD)) which are com¬ 
monly available for high-redshift galaxy samples, and 
E{B — U)gas, which can be difficult to constrain in 
the abs ence of Balmer decrement measurements. Fig¬ 
ure shows the difference between gas and contin¬ 
uum color excess {E{B — U)gas — E{B — U)stars), and 
the difference between the attenuation of recombination 
and continuum photons (Aho, — Acont): as a function 
of the SFD-determined sSFR, sSFR(SED). Note that 
sSFR(SED) combines the SFD-determined SFR with the 
SFD-determined stellar mass, so it is not the same as the 
default sSFR assumed throughout the paper (i.e., the 
latter combines SFR(Ha) with M*). We chose to quan¬ 
tify these relationships in terms of sSFR, and not SFR, 
as the correlations of E{B — U)gas — E{B — U)stars and 
Aho — Acont with sSFR(SED) are stronger than those 
with SFR(SED). To establish useful empirical relations 
between these quantities, we fit them with functions of 
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SFR(SED) (Mg yr-i) - MOSDEF 



Figure 22. Left: Comparison of SFR(Ha) assuming the |Cardelli et al. ||1989|l extinction curve (and the Balmer decrement), and SED- 
based SFRs assuming the MOSDEF attenuation curve. Higfit (bottom): Hatio of the dust-corrected Ho and UV-based SFRs, plotted 
against SFR(Hq:). The UV-based SFR is calculated from the flux of the best-fit SED model at 1600 A, and corrected for dust assuming 
E{B — V)stars- The red line and shaded region indicate the best-fit third-order polynomial fit and 68 % confidence interval between 
SFR{Hq;)/SFR(UV) and log(SFR(HQ;)/M 0 yr~^). The large stars in both panels denote the average values for the H/3-undetected galaxies. 
Right (top): Fraction of galaxies where SFR(Hq;) exceeds the UV-based SFR by more than Icr, in bins of SFR(Ha). 



Figure 23. Ratio of the observed Ho; and UV SFRs (i.e., uncor¬ 
rected for dust attenuation) versus the dust-corrected Ha SFR. The 
large star denotes the average values for the H/3-undetected galax¬ 
ies. Those galaxies with dust-corrected SFR(Ha)< lOMQyr”^ 
generally have observed SFR ratios that are consistent with unity 
within 3 cr. 

^ = l./(log[sSFR(SED)/yr-i] + 10): 
E{B-V)^^-E{B-V) 

stars — 0.049+ 0.079/C, (13) 

and 

Ana - Aeont = -0.0035 + 0.1972/e. (14) 



sSFR(SED) (yr-i) 


The scatter in these relations are a[E{B — V)gas — E{B — 
l^)stars] « 0.12 and cr[AHa “ Acont] « 0.24. Given 
E{B — V)stars and sSFR(SED) obtained from standard 
SED-fitting, one can use Equations and to then 
roughly estimate E{B — V)gas and Ana, keeping in mind 
the large scatter in these quantities. Einally, note that 
the difference in color excesses is found to decre ase on 
age wi th increasing sSFR(SED) (e.g., see also Price 
.|2014|) , opposite of the trend that we found between 


averae 


et al. 


B{B — V )gas — E{B — V)stars and the default sSFRs that 
assume Ha-based SERs (Eigures [I^ and [T^. This be¬ 
havior simply reflects the systematic difference between 
SED-based and Ha-based SERs, as discussed above. 


Figure 24. Difference between the gas and continuum color ex¬ 
cesses (top), and difference between the attenuation (in magni¬ 
tudes) of recombination line and continuum photons (bottom) ver¬ 
sus sSFR(SED). Here, sSFR(SED) is computed using the SED- 
based SFR and M*. The solid lines and shaded regions denote the 
best-fit relations and the 1 cr scatter, respectively, as indicated in 
the text. 


4.4.3. Summary of SFR Comparisons 

The physical picture put forth in the previous section 
leads to the expectation that UV- (or SED-) based SFRs 
should diverge from Ha-based ones—which we consider 
to be more reliable than the former—for galaxies with 
large SFRs, as an increasing fraction of the nebular line 
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and bolometric luminosities arise from regions of higher 
dust covering fraction. The statistics indicating an in¬ 
creasing SFR(Ha)/SFR(UV) ratio with SFR(Ha) are 
consistent with this scheme. A similar effect is seen in 
the comparison of UV and IR-based SFRs (e.g.JG oldader 
et al.|2002t |Reddy et al.|2010||Howell et al.||20lup, where 
the former generally underpredict the latter for galaxies 
where most of the star formation occurs in optically-thick 
regions that are essentially “invisible” in the UV. Be¬ 
cause the attenuation curve derived for the continuum is 
most sensitive to the less-reddened stellar populations in 
a galaxy, it is not surprising that inferences of total SFRs 
based on the UV continuum emission will underestimate 
the total SFR for very dusty galaxies. The assumption 
of UV- or SED-based SFRs may also lead to erroneous 
conclusions regarding the dependence of the difference in 
color excess between the stars and ionized gas, and the 
sSFR (e.g., Figures fl^ and [^. 

Thus, the common practice of assuming either equal 
color excesses of the ionized gas and stellar continuum, 
or assuming that the ionized gas color excess is sys¬ 
tematically larger than that of the stellar continuum, 
may lead to incorrect inferences of the SFRs. We also 
highlight the importance of treating the nebular emis¬ 
sion and stellar continuum with the appropriate atten¬ 
uation/extinction curves. Here, we have made indepen¬ 
dent measurements of the color excess of the ionized gas 
(aided by the Balmer decrement), and have calculated 
the the attenuation curve relevant for the stellar con¬ 
tinuum. Based on these analyses, we find disparities be¬ 
tween the Ha- and UV-based SFRs that signify a varying 
contribution of dusty stellar populations to the bolomet¬ 
ric luminosities of high-redshift galaxies. 

5. CONCLUSIONS 

We use a sample of 224 star-forming galaxies at red- 
shifts z = 1.36 — 2.59 with measurement of Ha and 
H/3 emission lines from the MOSDEF survey to inves¬ 
tigate the effects of dust on the stellar continuum of 
high-redshift galaxies. Our analysis suggests that the 
reddening of the UV continuum with increasing Balmer 
line opacity is most directly tied to dust obscuration. We 
calculate the attenuation curve by constructing, and tak¬ 
ing ratios of, the composite photometric measurements of 
galaxies in bins of specific SFR and Balmer optical depth. 
The attenuation curve derived here is very similar in 
shape to that of the SMC curve at A > 2500 A. Blueward 
of this wavelength, the curve has a shape that is identical 
to, but with a n ormalization that is A fc(A) ~ 1.63 lower 
than, that of the Calzetti et al. (2000) attenuation curve. 
Moreover, we find a marginal trend between the depth of 
absorption at 2175 A and the color excess of the stellar 
continuum. 

We explore the implications for the attenuation curve 
on the stellar populations and SFRs of high-redshift 
galaxies. In general, the new attenuation curve implies 
slightly redder color excesses {E{B — V)), « 20% lower 
SFRs, and stellar masses that are A\og{M*/M q) « 
0.16 dex lower, than those computed us ing the typi¬ 
cally assumed starburst attenuation curve (Calzetti et al. 
2000). We also find that while the color excess of the 


ionized gas is similar to that of the stellar continuum for 
about half of the objects in our sample, the difference 
in the reddening—and in the total attenuation—of the 


gas and stars in high-redshift galaxies is a strong func¬ 
tion of SFR. In particular, we find that the ionized gas 
is more reddened relative to the stellar continuum with 
increasing SFR. Our results are consistent with a simple 
picture in which the UV through optical continuum of 
high-redshift galaxies is dominated by a modestly red¬ 
dened stellar population, while a second, dustier, stellar 
population begins to dominate the bolometric luminos¬ 
ity and nebular line luminosity with increasing SFR. In 
this picture, UV- and SED-based SFRs may diverge from 
(i.e., underpredict) the total SFR even for galaxies with 
relatively modest Ha-based SFRs (> 2OM0yr“^). 

Our analysis and results demonstrate the utility of us¬ 
ing high quality near-IR spectroscopic data, along with 
multi-wavelength photometry, to directly measure the at¬ 
tenuation curve at high redshift. Further spectroscopy 
from the ongoing MOSDEF survey, along with deep mid- 
and far-IR imaging data, will allow us to build upon these 
results, and will be crucial for removing a key uncertainty 
(namely the shape and normalization of the attenuation 
curve) in the measurement of star-formation rates and 
the interpretation of stellar populations in the distant 
Universe. 
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